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I RANDOM  WATER  WAVES  UN HE ACMES 

* C.  Y.  Yang  and  Y.  H,  Chen* 


1 . Introduction 

This  paper  p^'esents  the  mathematical  derivation  and  physical  interpretation 
of  the  formula  governing  the  change  of  energy  as  random  waves  propagate  on  beaches 
from  deep  water  towards  the  shoreline.  The  energy  change  is  characterized  by  the 
wave  spectral  density.  The  problem  is  limited  to  linear,  one-dimensional  steady 
waves  directed  normally  to  a beach  of  constant  slope.  In  the  wave  field  it  is 
furthur  assumed  that  there  is  no  wave  reflection  from  the  shoreline;  no  energy  input 
from  wind  and  dissipation  from  turbulence  and  no  wave-wave  interaction. 

In  five  parts  the  paper  first  reviews  some  relevant  results  of  the  deter- 
ministic wave  theory  and  then  examines  the  classical  result  of  Longuet-Higgins  (1) 
on  the  transformation  of  a continuous  wave  spectrum  based  on  the  concept  of 
the  conservation  of  energy  flux  and  the  concept  of  random  noise  by  Rice  (2) . 

To  improve  on  Longuet-Higgins'  significant  but  approximate  result  where  the 
local  slope  of  the  beach  is  assumed  to  be  horizontal,  the  more  accurate  theory 
by  Friedrichs  (3)  for  deterministic  waves  is  used  to,  study  the  random  wave  problem. 

The  basic  idea  is  to  determine  wave  energy  changes  directly  from  the  exact 
solution  for  the  wave  surface  p.  This  simple  approach  is  possible  because  the 
exact  solution  n for  the  sloping  beach  gives  the  spatial  variation  in  wave  ampli- 

i 

tude  which  is  not  provided  by  the  wave  solution  for  a flat  bottom.  It  is  shown  i 

that  for  small  beach  slopes,  the  asymptotic  solution  of  Friedrichs  leads  to 
the  formula  governing  the  transformation  of  random  wave  spectrum  that  is  in 
complete  agreement  with  that  of  the  approximate  theory  by  Longuet-Higgins. 

For  beach  slopes  that  are  not  small,  a new  formula  for  the  wave  spectrum 
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transformation  is  derived  by  extending  Stoker's  (4)  exact  solution  for 
deterministic  waves  to  a random  wave  field.  Numerical  results  based  on  the  new 
formula  are  obtained  and  presented  along  with  those  by  Longuet-Hl ggins . 

2 . Review  of  Deterministic.  Water  Waves 

Consider  a one-dimensional,  linear,  progressive  single  wave  in  water  of 
uniform  depth  h.  The  wave  surface  with  amplitude  a,  wave  number  k,  frequency 
o,  and  phase  angle  0 may  be  denoted  by 

ri(x,t)  = a Cos(kx  + at  + 0)  (1) 

The  wave  energy  per  unit  horizontal  surface  area  Is 

E = Y pga^  (2) 

When  such  a single  wave  steadily  propagates  from  deep  water,  in  a normal 
direction,  towards  the  beach  with  a uniform  beach  slope  and  when  it  is  assumed 
that  there  is  no  energy  input  from  wind  and  no  dissipation  due  to  turbulence 
within  the  wave  field,  then  the  energy  flux  must  be  conserved.  That  is, 

EC  = E C (3) 

00  00 

where  C is  the  group  velocity  and  the  subscript  °°  denotes  the  quantity  in  deep 
water. 

The  above  results  for  a single  wave  can  be  extended  to  a multiple  of  n 
single  waves.  Let  the  combined  wave  surface  be  denoted  by  c(x,t).  Then  the 
wave  surface  and  the  energy  equations  are 

i;(x,t)  = E a Cos(k  X + 0 t + 0 ) (4) 

n nnn 
n 

E = I pg  E a^^  (5) 
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wJiere  the  subs^'ript  n denotes  the  nth  component  wave. 

3 . LonRuet-HIgBj-ns' Approximate  Theory  (1957) 

Consider  a random  progressive  wave  surface  n(x,t)  denoted  by  the  same 
equation  (1)  except  that  in  this  case  the  phase  angle  6 is  a random  variable, 
distributed  uniformly  in  (0,2ir).  Thus  n(x,t)  is  a random  process.  Note  that 
for  this  random  wave  model,  however,  Eqs.  (2)  and  (3)  remain  valid  since  the 
energy  and  the  group  velocity  for  each  wave  sample  are  Independent  of  the  random 
phase.  All  wave  samples  have  equal  energy  E and  group  velocity  C so  that  these 
quantities  can  be  considered  as  ensemble  mean  values.  The  ensemble  mean  for  the 
random  process  n(x,t)  is  zero  and  the  variance 
Var[n(x,t)]  = n^(x,t)  = Y 

where  the  bar  over  the  symbol  denotes  ensemble  average.  Consequently,  in  view 
of  Eqs.  (2)  and  (6)  the  variance  or  the  ensemble  mean  square  surface  represents 
the  wave  energy  E except  for  the  constant  pg.  This  connection  leads  to  the 
familiar  concept  of  an  energy  spectrum  density.  Thus  if  S(k)  is  defined  by 

S(k)Ak  = Var[n(x,t)]  - E (7) 


it  represents  the  energy  density  associated  with  random  waves  of  wave  number  k, 
simply  the  wave  number  spectrum.  Using  the  equation  for  the  conservation  of 
energy  flux,  Eq.  (3),  yields  the  following  formula, 

S(k)AkC  = S (k  )Ak  C (8) 

00  00  00  oo  ' ' 


If  it  is  assumed  that  the  frequency  o is  constant,  then  the  ratio  of  the  group 
velocities  C to  C is  equal  to  that  of  the  wave  number  increment  Ak  to  Ak 

CO  OO 
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and  In  this  way  hongueC-Higy.ins ' important  formula  may  lie  obtained  (1)  as 


S(k)  = S^(kJ  (9) 

for  the  conservation  of  wave  number  spectrum.  This  conservation  formula 
for  the  single  random  wave  model,  Eq.  (1),  with  a single  random  pliase  angle, 
0,  applies  to  a multiple  of  n random  waves  defined  by  Eqs.  (4)  and  (5), 
which  can  be  shown  as  follows. 

Let  it  be  assumed  that  the  random  phase  angles  6^  in  Eq.  (4)  are 
Independently  and  uniformly  distributed  in  (0,27r).  Then  for  the  combined 
wave  surface  the  ensemble  mean  is  zero  and  the  variance 

Var[i;(x,t)]  = ?^(x,t)  = (10) 

l n 
n 

so  that  the  wave  energy  associated  with  the  component  random  wave  of  wave 
number  is  S(k^)Ak.  It  is  again  governed  by  the  relation  for  the  conser- 
vation of  energy  flux,  Eq.  (8),  but  under  the  additional  assumption  of 
independent  propagation  for  all  component  waves.  Finally,  under  the  same 
assumption  of  constant  frequency  for  each  component  wave,  the  same  formula 
Eq . (9)  for  the  conservation  of  the  wave  number  spectrum  holds. 

Thus,  It  is  shown  that  under  the  assumption  of  no  wave-wave  interaction, 
the  spectrum  transformation  relation  for  the  random  multiple  wave  model  can 
be  derived  from  a single  random  wave  model.  This  simplification  is  used  in 
the  following  analysis. 

The  theory  is  considered  approximate  because  the  basic  wave  surface 
Eq.  (1)  and  the  group  velocity  C are  derived  under  the  assumption  of  a flat 
bottom.  This  assumption  is  removed  in  the  following. 
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4.  Random  Waves  on  a Shallow  SlopinR  Beach 

For  a shallow  sloping  beach,  Friedrichs'  asymptotic  solution  (3,  p.  114)  for  a 
progressive  wave  is  given  by,  (Figure  1), 

n(x,t)  = 2k  ^ A(X)Sin{a)  ^K(X)  + at  + v)  (11) 

00  oo  ^ 

where  w is  the  small  beach  angle,  X is  the  ratio  of  pseudo  wave  length  at  x 
to  that  at  X -»■  «>  and 


uxk  = X tanh  ^ X 


A(X)  = 


I 2ttu) 


■ii 


(X) 


J(X)  = X + (1  - X^)tanh  ^ X 


K(X)  = tanh  ^ X + 


, -1  dv 
tanh  V — 

V 


In  deep  water,  wx  ->  «> 


(12) 

(13) 

(14) 

(15) 


“ 1 — ^ 

Ti\X,t)  = 2k  72710)  Sin(k  x + ot+^  -t)  (16) 

The  term  pseudo  wave  length  is  used  because  the  wave  solution  Eq.  (11)  does  not 
have  a wave  number  term.  Consequently  only  an  approximate  or  pseudo  wave 
length  can  be  defined.  Using  such  a definition  for  X in  Eq.  (12)  leads  to 
the  well-known  dispersion  relation  for  small  beach  angle,  ui. 

Since  the  governing  equation  for  the  wave  field  and  the  boundary  conditions  are 
linear  and  homogeneous,  a phase  angle  P can  be  introduced  in  Eqs.  (11)  and  (16) 
to  form  the  following  modified  solution. 


= 2k  ^ A(X)  Sin{(i)“^  K(X)  + a t + |-  + 6} 


n(x,t) 


(17) 


and  for  lox  -* 


n(x,t)  = 2k.  ^ Jl-nui  Sin{k  x + a t + ~ 

O)  OOJ 


- t - 0) 


(18) 


Consider  now  that  the  phase  angle  is  random  with  uni  form  distribution  in 
(0,  2n) . Ihen  the  variance,  the  spectruin  and  the  transformation  for  the  spectrum 
can  be  derived  as  follows.  From  F.q . (17)  and  by  definition, 

Vartn(x,t)l  = 4k“^  A^(A)  = S(k)Ak  (19) 

From  Eq . (18) 


Var[n(x,t)]  = 8k^^TTu  = S^(k^)Ak^ 


(20) 


Dividing  Eq.  (19)  by  Eq.  (20)  yields 
S(k)  = S (k  ) ~ 

oo'  CO  2710)  dk 

Using  the  definitions  of  A,  A(X) , j(A)  and  Eq.  (12),  the  above  equation  is 
reduced  to 


(21) 


S(k)  = S„(k J 


(22) 


Thus  it  is  shown  that  based  on  Friedrichs'  asymptotic  solution  for  a shallow 
sloping  beach,  the  transformation  relation  for  wave  number  spectrum  can  be  derived 
directly  without  the  assumption  of  a flat  bottom.  The  result  (Eq.  (22))  for 
shallow  slopes  agree  completely  with  that  from  Longuet-Higgins ' approximate 
theory,  Eq.  (9).  The  condition  of  the  energy  flux  conservation  and  thus  the 
group  velocity  are  not  needed  in  the  derivation. 

5 . Random  Waves  on  Non-Shallow  Sloping  Beach 

For  a non-shallow  sloping  beach  of  45°,  Stoker's  exact  solution  (4,  p.  24) 
are  given  by  two  potential  functions,  a regular  function  i(7j(x,0,t)  and  a singular 
function  i(/2(x,0,t), 
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II 


+ Cos  X 


Sin  X J 


(23) 


It  , -X 

e [e 


It 


i|i2(x.0,t)  - {C  (x)  [Sin  x - Cosx] 

Fl  ^ 


[-2  + S^(x)] 


[Cosx  + Sin  x]  - e ^E.(x)) 


(24) 


where  for  simplicity  all  quantities  in  the  above  two  equations  and  the  following 
equations  up  to  Eq.  (30)  are  dimensionless  with  x defined  as  the  product  of 
the  horizontal  dimensional  coordinate  and  the  wave  number,  k : and  t as  that 
of  the  dimensional  time  and  the  frequency  S^(x),  C^(x)  and  E^(x)  are  Sine, 
Cosine,  and  exponential  integral  functions,  respectively. 

From  the  potential  functions  two  standing  wave  surface  solutions  are 

selected.  These  are: 

9iii 


n^(x,t)  = = -i|i^(x,0,0)Sin  t 

_ _7T 

3i|'2 

n2(x,t)  = = il'2(x,0,0)Cos  t 

1 


(25) 


7T  “X 

— (e  + Cos  X - Sin  x)Sin  t 


(26) 


{C.  (x)  [sin  X --  Cos  x]-[y  + S.  (x)]  [Cos  x + Sin  x] 


-e  E^(x) }Cost 


These  standing  waves  can  be  combined  as  follows,  using  ijij^(x)  = iji^(x,0,0)  and 

iji2(x,0,0)  = ij;2(x)  , 

n(x,t)  ^ [—i  (x)Sint]  - - [iii„(x)Cost] 

' H 1 IT  Z 


(27) 


which,  as  X ^ reduces  to  a simple  progressive  wave  given  by 
n(x.t)  = — (Cos  X - Sin  x)Sint  + — (Cos  x + Sin  x)Cost 

n /2 


X-w 


(28) 


Sin(x  + t + ■2I') 
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As  in  the  previous  section,  a phase  angle  0 is  now  introduced  to  obtain 


the  modified  wave  solutions, 


n(x,t)  = ii)^(x)Sin(t  + 0)  - ^ ijj,(x)Cos(t  + 0)  (29) 

oo  " Sin(x  + t + |-  + 0)  (30) 

Having  established  the  wave  solutions  in  proper  forms,  the  phase  angle  0 
is  now  considered  to  be  a random  variable  with  uniform  distribution  in  (0,  2n) . 
Then  from  Eq.  (29)  and  by  definition. 


Var[n(x,t)]  = ^ ^ li^^^x)]  - S(k)Ak 

it''  Tt 

From  Eq.  (30) 

Var[n(x,t)  ] = 4 = S (k  )Ak 

’ 2.  coco  oo 

Hence  S (k  ) B(x)  -7^ 

OO  00  QK. 

wher-  . 

TI 


(31) 

(32) 

(33) 

(34) 


The  transformation  relation  of  wave  number  spectrum  for  random  water  waves 
on  a non-shallow  sloping  beach  of  45°  is  thus  established  by  Eq . (33),  based  on 
Stoker's  exact  solution.  Again,  as  shown  previously,  the  single  random  wave 
solution  is  applicable  to  a sura  of  many  random  wave  components  under  the 
assumption  of  independent  phase  angles  0^  and  independent  v^ave  propagation. 

A comparison  of  Eq.  (33)  with  Longuet-Higgins ' solution,  Eq.  (9),  shows  that 
the  correction  to  the  approximate  theory  is  the  product  B (x) (dk^/dk) . 

Note  that  since  only  the  ratio  of  the  variances  is  relevant  in 
Eqs.  (31)  and  (32),  S(k),  Ak,  S^(k^) , Ak^  and  all  quantities  to  be  used 
hereafter,  with  the  exception  of  t}ij(x)  and  (j;2(x),  can  be  considered  as 
dimensional  for  simple  physical  interpretation. 
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6.  Numerical  Results 


Since  the  dispersion  relation,  = k tanh  kh  where  h is  the  water  depth, 
is  no  longer  applicable  in  the  exact  solution,  it  is  very  difficult  to  calculate 
accurately  the  term  dk^/dk  in  Eq.(33).  Furthermore,  in  practical  applications, 
the  wave  frequency  spectrum  S(ri)  is  generally  preferable  over  the  wave  number 
spectrum.  Therefore,  it  is  useful  to  examine  the  alternate  form  S(a). 

Using  the  frequency  spectrum  in  Eqs.  (31)  and  (32)  yields 

Var[ri(x,t)]  = j [i|j^(x)  + ] = S(a)Aa  (35) 

2Tr 

Var[ri(“>,t)  ] ^ J = S^(a)Aa  (36) 

where  the  same  assumption  a = is  used  as  before.  From  these  equations,  the 
alternate  relation  for  wave  energy  transformation  is  obtained  as 

S(o)  = ~ [i|^^(x)  + <|j2(x)l  S^(a)  = B(x)S^(o)  (37) 

IT 

This  form  is  clearly  superior  than  that  with  the  wave  number  spectrum  S(k), 
for  the  troublesome  term  dk  /dk  is  no  longer  needed.  In  terms  of  S(a), 

oo 

Longuet-Higgins ' approximate  solution,  Eq.  (9),  becomes 

S(o)  = ^ S (a)  = ^ S (o)  (38) 

C » dk 

00 

where  the  term  dk/dk  can  be  easily  calculated  according  to  the  dispersion 
relation  of  the  linear  wave  theory  for  a flat  bottom,  which  is  consistent 
with  the  approximate  theory. 

The  new  transformation  factor  B(x)  in  Eq . (37)  is  plotted  in  Figure  1 
together  with  the  factor  dk/dk^  in  Eq.  (38)  which  is  equal  to  C^/C  and  also 
to  E/E  according  to  the  approximate  theory. 

OO 
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7.  Conclusion 


( 

i 


A new  formula  Eq . (37)  is  derived  for  the  change  of  wave  energy  as  random 
waves  from  deep  water  propagate  toward  tlie  beacli.  This  formula  is  considered 
to  be  an  improvement  to  the  well-knoii;n  one  by  l,onguel-Hi  ggins  becatise  the 
assumption  of  a flat  bottom  locally  is  removed.  The  basic  idea  is  to  use 
the  spatial  variation  of  the  wave  amplitude  from  the  exact  solution  of 
Stoker  in  the  derivation  of  wave  energy  rather  than  the  ccjncept  of  con- 
servation of  energy  flux. 

For  a shallow  sloping  beach,  it  is  shown  that  the  formula  derived  from 
Friedrichs'  asymptotic  solution  agrees  completely  with  that  by  Longuet-Hlggins . 
For  a non-shallow  sloping  beach  of  45°,  numerical  results  are  computed  from 
the  explicit  new  formula  and  compared  with  that  by  Longuet-Higgins . The 
comparison  shows  a general  agreement  in  their  spatial  variations.  The 
quantitative  difference  between  them  is  within  about  8 percent. 
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LINEAR  RANDO^I  WAVES  ON  WATER  OF  NONUNIFORM  DEPTH 

SUMMARY 

Propagation  of  random  gravity  waves  on  water  of  variable  depth  Is 
examined  to  Incorporate  nonuniform  bottom  effects  explicitly  into  known 
random  wave  models  assuming  a locally  flat  bottom.  The  methodology  is 
the  same  WKB  approximation  previously  employed  by  Mei,  et  al.  [1968]  and 
Chu  and  Mei  [1970]  in  the  case  of  monochromatic  waves.  However,  attention 
is  restricted  to  the  linear  case  and  concepts  are  extended  to  incoherent 
random  waves  which  admit  an  orthogonal  spectral  representation.  This  is 
shown  to  be  valid  in  a slowly  varying  inhomogeneous  medium  just  as  in  a 
homogeneous  case.  Nonuniformities  in  depth  introduce  a skewness  in  the 
free  surface  whereas,  to  the  same  order,  the  corresponding  spectral  density 
remains  unchanged.  Refractive  transformation  of  continuous  spectra  is 
re-examined  emphasizing  a simplified  approach  for  applications.  This  is 
demonstrated  with  an  approximate  closed  form  solution  valid  for  spectral 
predictions  in  shallow  water.  Finally,  an  inhomogeneous  Gaussian  wave 
model  is  considered.  Associated  envelope  and  wave  height  distributions  are 
derived. 


I 

i 


, INTRODUCTION 

As  waves  propagate  from  a homogeneous  region  such  as  deep  water  into 
a region  with  uneven  underwater  topography,  their  characteristics  vary 
spatially.  Prediction  of  such  variations  has  been  systematically  investi- 
gated in  the  case  of  monochromatic  waves  by  several  investigators  using 
an  approach  known  as  the  WKB  approximation  [see,  e.g.,  Mei,  et  al.  1968; 

Chu  and  Mei,  1970].  The  basis  of  this  approach  is  an  asymptotic  expansion 
of  the  governing  equations  of  wave  motion  in  powers  of  a small  parameter 
e which  characterizes  the  slow  variation  of  the  still  water  depth.  Formally, 

1 

e is  proportional  to  the  fractional  changes  in  water  depth  over  a given 
wavelength  for  steady  state.  Partial  reflections  that  always  occur  at 
changes  in  the  medium  are  neglected.  However,  these  are  known  to  be 
exponentially  small.  This  approach  leads  to  the  derivation  of  a multitude 
[ of  results  on  the  spatial  variation  of  kinematic-dynamic  wave  properties 

in  a systematic  manner  directly  through  the  governing  equations  of  wave 
motion  in  contrast  with  previous  wave  theories  which  consider  the  bottom 
to  be  horizontal  locally,  dealing  with  the  actual  depth  variation  after- 
wards indirectly  through  conservation  of  energy  flux. 

The  motivation  here  is  to  apply  the  approximation  to  incoherent 
random  waves.  To  establish  the  corresponding  results  in  this  case,  ideas 
due  to  Mel,  et  al . [1968]  and  Chu  and  Mei  [1970]  will  be  freely  drawn  on. 
However,  attention  is  focused  to  linear  waves,  and  the  primary  emphasis 
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is  probabilistic  with  a view  to  establishing  a spatially  Inhomogeneous 
random  wave  model.  In  the  first  section,  various  definitions  and  formulation 
of  the  WKB  approximation  is  summarized,  closely  following  the  vjork  of 
Chu  and  Mei  [1970]  with  the  free  surface  and  the  appropriate  velocity 
potential  expressed  in  a continuous  spectral  form.  Subsequently,  spatial 
variations  of  random  amplitudes  and  free  surface  spectra,  and  the  explicit 
effect  of  an  uneven  bottom  on  the  free  surface  and  its  probability  struc- 
ture are  examined.  Finally,  the  probability  distribution  of  wave  heights 
is  derived  for  a Gaussian  free  surface  without  any  Inherent  bandwidth 
constraints  on  the  associated  spectral  density. 

DEFINITIONS  AND  ANALYSIS 

Consider  the  steady  state  linear  wave  propagation-refraction  problem, 
using  a coordinate  system  in  which  x and  y represent  t}ie  horizontal  cartesian 
coordinates  fixed  at  still  water  level  with  the  z-axls  pointing  upwards. 
Letting  <|)(x,z,t),  Ti(x,t)  and  h(x),  where  ^ = (x,y),  denote  respectively  the 
velocity  potential,  the  free  surface  and  the  still  water  depth,  the  problem 
is  governed  by 

v2())  + = 0 , -h(x)  _<  z j<  0 

gn  + <!>(.  = 0 , z = 0 

nj.  + <(>^  = 0 , z = o 

4>  + Vh  • 
z — 


V(|)  = 0 


z = -h(x) 


(1  a,b,c,d) 
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In  the  above,  2.  ~ 3/3y,  0)  represents  the  horizontal  gradient  operator, 

and  ^ It  is  assumed  that  the  still  water  depth  h(x)  varies  slowly 

over  the  distance  of  a given  wavelength,  X.  This  spatial  variation  is 
formally  characterized  by  a parameter  c such  that 

c *=  0(l^h|/kh)  <<  1 (2) 

where  k = 2tt/X  is  the  wave  number.  In  this  manner,  we  can  Introduce  the 
following  slowly  varying  variables: 

X = Cx,y)  = ex  , z = z , t = et  (3) 

The  governing  equations  can  now  be  rewritten,  after  combining  (lb)  and  (Ic) , 
in  the  form 

+ <{,-_=  0 , -h(x)  1 z £ 0 

'*’tt  ” ~ ® ’ z = 0 

4,.  + e^Vh  • ^(t>  = 0 , z = -h(x)  (A  a,b,c) 

where  Y.  ~ “ (3/3x,  3/Dy,  0). 

A random  wave  field  which  is  stationary  in  time,  t,  and  inhomogeneous 
in  the  horizontal  space,  3c,  can  be  represented  by 

n(x, t)  = / dA(k,u),x)  exp(ix/e)  (5a) 

jc,a) 


with  an  associated  velocity  potential  in  the  form 
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<}i(x,z,t)  = / dB(k,w,x,z)  exp(lx/E) 

k.,a) 


(5b) 


The  phase  function 


X(x 


r- 

,t)=/  Ji'djc- 


wt 


(6) 


defines  the  horizontal  vector  wave  number  Jc  = (k2^,k2)  and  frequency,  : 


k(x)  = ^ , w = -x~ 


(7  a.b) 


such  that 


Vxk  = 0 and  Vu)  = 0 


(8  a.b) 


It  is  noted  that  A(]c,a),2c^)  and  B(_k,u),5c,  z)  in  (5  a.b)  are  random  amplitude 
processes  which  are  of  central  importance  to  the  specification  of  the  wave 
field.  For  a spatially  homogeneous  wave  field,  both  A and  B become 
independent  of  x and  constitute  the  usual  Fourier-Stielj es  representations 
for  the  free  surface  q and  the  velocity  potential  (fi.  In  particular,  the 
covariance  matrices  of  the  increments  dA  and  dB  are  diagonal  entirely  as 
a consequence  of  the  homogeneity  condition.  This  so-called  orthogonality 
property  would  not  strictly  be  true  for  a spatially  inhomogeneous  wave 
field  unless  the  amplitudes  A and  B can  be  regarded  as  slowly  varying 
functions  of  x such  that,  corresponding  to  a gradual  variation  in  depth,  h. 


r 


the  fractional  changes  in  A and  B over  a given  wavelength  are  much  smaller 
than  unity.  The  orthogonality  property  is  then  approximately  valid  so 


E{dA(k,u,xl  dA*  Oc'.w'.x)}  = 5k  doj  , (k  = k',  u = ta') 


otherwise  (9) 


where  £{•}  denotes  a probability  or  ensemble  average  for  the  argument,  and 
6_k  = dkjdk2  as  a shorthand  notation.  The  function  iJj  is  the  inhomogeneous 
spectral  density  representing  locally  the  density  of  contributions  to  the 
mean  square  surface  deformation,  l.e.. 


E{ln|2}  = / 6_kd 


It  is  assumed  that  the  random  differentials  dA  and  dB  admit  the  following 
expansions  of  WKB  type  [see,  e.g.,  Chu  and  Mei,  1970]: 


dA  = \ dA^^^(k,w,x)  , dB  = \ dB^^\k,u,x)  (11  a,b) 

j=0  j=0 

where,  for  a fixed  j,  dA^^  and  dB^^^  are  related  to  one  another  through  (lb) 
in  the  form 


dA^^^  (k,a),x)  = i(a)/g)  dB^^\k,a),x,0) 


(12) 
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Substitutions  of  (11)  and  (5)  into  (4  a,b,c)  yields  a series  of  pertur- 
bation equations  for  the  random  differentials  For  j = 0 and  1,  to 

which  the  analysis  will  be  restricted,  we  have 


zz 


-h  < 2 < 0 


gdBp^  - w^dB^J^  = 0 


2 = 0 


dB^^^  = 


(13  a,b,c) 


where 


= 0 


R^^^  = -1  { 2k  • VdB^^^  + V • k } 


= _1  Vh  • { k dB^°^  }_ 

— — 2 “ -h 


(lA  a,b,c) 


The  corresponding  solutions  to  (13a)  and  (13c)  are  In  the  form 


ni 

dB''-J''(k,io,x,2)  = dC^-^''  cosh  Q + ^ sinh  Q 


y'*Q 

R^^^  cosh  Q'  dQ'’ 

0 


- cosh  Q / R'^''  slnh  Q"  dQ""} 

0 
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where  Q = kC^  + h)  and  k = |jc| . On  substitution  from  (15),  (13b)  becomes 


(k  sinh  q - k cosh  q){dC 


+ (k  cosh  q - sinh  q)  {-j-  + 

K 


sinh  Q dQ) 


cosh  Q dQ} 


in  which  k^  = w^/g  and  q = kh.  In  essence,  (16)  provides  the  explicit  form 
of  the  differential  coefficient  dC^\  For  j = 0,  the  system  (13  a,b,c)  is 
homogeneous.  In  this  case  (lb)  leads  to  the  dispersion  relation 


k tanh  q = k = u^/g. 


which  requires  that  dA^^^  and  be  expressed  as 


dA^°^  = dA^^^(k,x) 


jn(0)  i^CO)  , „ ^ ,.(0)  cosh  Q 

dB  = dC  cosh  Q = -i(g/u)  dA  ; — ^ 

cosh  q 


(18  a,b) 


For  j = 1,  and  using  (14  b,c)  and  (15),  the  solution  of  the  inhomogeneous 
system  (13  a,b,c)  yields  dB^^^  in  the  form 


dB^^^  = dC^^^  cosh  Q - i dB^°^  F(k,x,5) 
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where 

F(k,x, z)  = tanh  Q + Q) 


with 


i' 

(20)  ; 


(21  a,b,c) 


In  this  case,  by  virtue  of  the  dispersion  relation  (17),  (16)  does  not 
provide  a unique  choice  for  the  differential  coefficient  dC^^\  This  coeffi- 
cient must,  therefore,  be  chosen  to  give  the  proper  limit  in  deep  water,  i.e., 
as  q ^ we  must  have  dB^^^  = 0 for  n ^ 1.  It  can  be  shown  then  [see, 
Appendix  A]  that  this  requirement  leads  to  the  conclusion  dC^^^  = 0,  and 
(19)  becomes 


dB^^^  = -i  dB^°^  F(k,x,S)  (22) 

The  preceding  solution  for  dB^^^  differs  from  the  corresponding  one  derived 
by  Chu  and  Mei  [1970],  It  is  believed  here  that  this  is  due  to  an  erroneous 
conlcuslon  on  their  part  as  to  the  limiting  values  of  the  quantities 
involved  in  (19)  as  q 

With  the  solution  for  dB^^^  at  hand,  it  is  Immediate  from  (12)  that 
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= -i  FQc.x.O) 


(23) 


and,  therefore,  up  to  the  order  O(e^): 


nCx,t)  = / 


[1  - ±E  F(k.,x,0)J  dA^*^^  exp  (ix/e) 


<|.(x,5,t)  = -ig  / [1  - ie  F(k,x,2)]  cos^q^'"^  expdx/O  (24  a,b) 

k 

It  will  be  expedient  from  now  on  to  work  with  the  original  variables  (x,z,t). 
Noting  that  ^ = e^,  we  define  = eS^  (j  = 1,2,3),  F(k,x,z)  = e F(k,x,z),  and 
(24  a,b)  become 


(0) 


n(x,t)  = / [1  - iF(k,x,0)]  dA^^^(k,x)  expdx) 
k 


,t)  = /[I  - 
k 

<^(x,z,t)  = -ig  /[I  - iFCk,x,z)]  exp(ix)  (25  a,b) 


where 


xCx.t) 


X 

= J k • dx 


- (Ot 


(26) 


Now,  by  using  a unit  vector  je^  in  the  z-direction,  other  relevant  properties 
of  the  wave  field  can  be  sunnnarized  up  to  0(|^h|^)  as  follows: 

Velocity  Field 

+ ii  e = f {VdB^°^  + ik[l  - iF(k,x.z)]  dB^°^  + dB^°^e3}  exp (lx) 


(27) 
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Acceleration  Field 


^ (^'}>  '*"‘^£3)  “ .f  + ljc[  1 - lF(k,x,z)]  dB^^^  + dB^*^^£2}exp(ix) 


(28) 


Free  Surface  Gradient 


^n(x, t)  = J' + ijc[l  - iF(_k,x,0)]  exp(ix) 


(29) 


It  is  evident  in  the  above  results  that  the  differential  coefficient  dA^^^ 
— and,  therefore,  dB^*^^—  is  so  far  arbitrary  in  terras  of  a probability 
structure  and  the  explicit  dependence  on  x-  In  the  following,  we  will  first 
discuss  the  spatial  dependency  of  dA^*^^  and  the  associated  free  surface 
spectra  since  this  can  be  done  in  a manner  independent  of  any  explicit 
probabilistic  constraints. 


SPATIAL  DEPENDENCY 

The  nature  of  x~dependency  of  the  random  amplitude  dA^^^  is  determined 
by  requiring  that  dB^*^^  and  dB^^^  satisfy  the  consistency  or  solvability 
condition  for  the  inhomogeneous  system  (13  a,b,c).  This  condition  corresponds 
to  (16),  which  becomes,  on  substitution  from  (14  b,c) 


Vh  • [k  dB^°^]  , + 

— — ■*  7.  =5  — h 


if! 


dB^°^  + dB^°^  V • k}  cosh  Q dQ  = 0 


(30) 
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Making  use  of  (18b)  and  carrying  out  the  integration,  we  obtain 


k 

k 


V dA 


(0) 


dA 


(0) 


2 (q  tanh  q - 1)  ^ 

slnh  2q  + 2q  k 


Vq 


(31) 


In  principle,  this  equation  together  with  the  irrotationality  condition  on 
k,  i.e.  (8a),  and  boundary  conditions  associated  with  k and  dA^^^  provides 
the  description  for  the  spatial  variitlon  of  dA^^\  However,  realizing 
that  one  is  in  fact  dealing  with  an  ersemble  of  amplitudes  dA^*^\  a numerical 
approach  which  is  Implicitly  required  by  this  description  in  its  present 
form  above  is  highly  undesirable.  At  the  same  time,  a careful  examination 
of  C31)  suggests  that  it  can  be  transformed  further  by  a change  of  variables 
from  X = Cx,y)  to  a curvilinear  system  (s,n)  where  s and  n denote, 
respectively,  coordinates  measured  along  the  orthogonal  trajectories  (rays) 
and  the  crests  of  a component  wavelet  as  schematically  illustrated  in 
Figure  1.  This  transformation  Isee,  Appendix  BJ  simplifies  C31)  to  a form 
which  can  be  integrated  to  give 


dA«> 


g 


(k,x)  = constant 


C32) 


along  a wave  ray.  In  the  above,  p represents  the  separation  distance  between 
adjacent  rays,  and  C is  the  group  speed  given  in  general  by 

O 


c 

g 


9u)  _ _w  i slnh  2q  + 2q 
9k  2k  ( sinh  2q 


(33) 
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It  is  evident  that,  for  waves  propagating  from  a homogeneous  region  such 
as  deep  water  into  a spatially  inhomogeneous  regime,  (32)  implies 


(C  )l'2dA(«(k  ) 

*^co  ® 


pl/2  ^in 


where  the  subscript  (“)  designates  the  constant  values  in  the  homogeneous 
region.  It  is  immediate  from  (34)  that  we  have  the  desired  relation 


(0)  ( i (0) 

dA^\k,x)={-}  dAl^^k^) 


to  describe  the  spatial  modulation  of  the  amplitude  dA^^\  Obviously  (35) 
suggests,  first,  that  dA^*^^  is  linearly  proportional  to  dA^^\  However, 
since  the  proportionality  constant  is  real-valued,  both  dA^^^  and  dA^^^ 
have  the  same  phase.  Further,  since  the  same  constant  is  also  deterministic, 
the  probability  structure  of  dA^^^  is  readily  obtained  once  that  of  dA^^^ 
is  specified.  Finally,  realizing  that  (^)  is  a random  process  with 

orthogonal  increments,  the  most  fundamental  result  embedded  in  (35)  is  the 
proof  to  the  heuristic  argument  that  the  inhomogeneous  amplitude  A^*^^  and, 
therefore,  A^^\  etc.,  have  orthogonal  increments  if  the  medium  is  slowly 
varying. 

At  this  point,  we  may  note  that  the  ratios 


(36  a,b) 
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are  recognized  as  the  ray  separation  factor  and  the  refraction  factor, 
respectively.  The  spatial  variation  of  0 implicitly  required  in  (35)  is 
governed  by  the  equation  of  wave  intensity,  wliich  has  been  previously 
discussed  in  detail  by  Munk  and  Arthur  [1952],  and  need  not  be  further 
elaborated  on  here.  It  suffices  to  point  out  that  the  determination  of 
0 in  general  requires  numerical  calculations,  although  approximate  solutions 
of  the  WKB  type  that  may  very  well  be  of  relevance  here  are  indicated 
also  by  the  same  authors. 


IIJHOMOGENEOUS  TREE  SURFACE  SPECTRA 


By  virtue  of  (9),  (11a)  and  (17),  the  inhomogeneous  free  surface 
spectrum  can  be  expressed  in  the  form 

^(k,x)  = i^^°^(k,x)  + i|^^^^(k,x)  + 0(|Vh|2)  (37) 

where 

i|j^*^\k,x)  = E{  |dA^°^  |2}/6k 

i/i^^\k,x)  = E{dA^°^  dA*^^^  + dA*^°^  dA^^h/5k  (38  a,b) 

Now,  using  (18b),  (23)  and  recognizing  that  F(k,x, 0)  is  real,  one  can 
verify  that 

«(;^^^(k,x)  = 0 (39) 


and,  (27)  becomes 
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1 

f 

f 

i 

i 


1. 


ij^Ck.x) 


<P^^\k,x) 


0(|Vh|2) 


(40) 


Therefore,  entirely  as  a consequence  of  the  fact  that  the  solutions  of  the 
two  leading  orders  dA^*^^  and  dA^^^  are  out  of  phase  by  (™/2),  the  free 
surface  spectral  density  is  characterized  by  corresponding  to  a 

locally  flat  bottom  correct  up  to  the  order  0(|_Vh|^).  Consequently,  the 
spatial  transformation  of  4;(_k,x),  which  is  of  main  interest  here,  is  described 
to  the  order  0C|^h|^)  in  terms  of  the  spatial  modulation  of  dA^^\  It 
follows  then,  by  multiplying  (32)  by  its  complex  conjugate  and  using  C38a) , 


P Cg  t|)^°\k,x) 


6k.  = constant 


(41) 


along  a wave  ray.  This  is  obviously  nothing  but  the  conservation  of  energy 
flux,  and  one  can  further  simplify  it  as  previously  done  by  Longuet-Hlgglns 
[1956]  to  obtain  the  classical  result  that 


4.<^°\k,x)  = /°\k  ) = 


constant 


(42) 


along  the  wave  ray.  Another  way  the  same  result  is  obtained  is  to  follow 

an  approach  that  corresponds  in  essence  to  what  Chu  and  Mel  [1970]  have 

done  in  the  case  of  monochromatic  waves.  In  this  case,  we  multiply  (CO) 

* 

by  (dA  '/cosh  q)  and  use  the  Leibniz  rule  to  obtain 

dB^®^  I 2 dz  --  0 (43) 

Making  use  of  (18b)  and  ensemble  averaging,  one  can  carry  out  the  Integra- 


I 


tion  to  obtain 
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V • [C  E{ I 2}]  = 0 


(44) 


where  C = k . Now,  using  (38a)  and  the  fact  that  V 
~8  S “ — 

e.g.,  Phillips,  1969,  p.  147],  (44)  becomes 


(C  6k)  = 0 [see, 
- 


k 


= 0 


(45) 


which  is  identical  to  (42),  recognizing  that  ji  * ^ = 3/3s. 

In  summarizing  the  preceding  discussion,  we  observe  that  the  spectral 
amplitude  remains  invariant  along  a wave  ray  although  the  associated 

wave  space  is  distorted  due  to  the  spatial  modulation  of  the  magnitude  and 
the  direction  of  Consequently,  the  complete  spectral  transformation 
requires,  in  addition  to  the  amplitude  information  a knowledge  of 

how  the  vector  wave  number  Jc  varies  along  a wave  ray.  This  is  embedded  in 
the  well-known  Hamiltonian  prope'ties  of  the  wave  field  [see,  e.g., 
Bretherton  and  Garrett,  1969],  which  can  be  expressed  for  the  steady  state 
case  of  interest  here  as 


dx 

3s 


and 


3k 

Is 


_1 

C 

g 


(Vm), 

— k = const 


where,  from  (17), 


(46  a,b) 


(Vio), 

— k = const 


ko) 


Vh 


sinh  2q  — 


(47) 
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It  is  seen  that  (46  a)  obviously  describes  a wave  ray,  and  (46a)  characterizes 
the  variation  of  along  the  same  ray.  In  essence,  (46  b)  represents  the 
kinematical  conservation  of  the  vector  wave  number  jc.  It  is,  therefore, 
a restatement  of  (8  a,b)  since  we  can  show,  by  recognizing  that  u = u(Jc,x), 
jk  = It (x) , and  expanding  (8  b)  with  ^ x k = 0,  that 


, 8k 

||._+(Vu.)k=0 


(48) 


where 


3a)  ^ = r 

8k  * 8x  “ 


3k 

8x 


~ 8k  ■ 

oCl)  — 

Ik 


w fv  - 

3y  _ 


Now,  noting  that  (8a)/8jc)  = C k and  Ic  • (3/8x)  = 8/8s,  (46  b)  follows. 

8 

Consider  now  random  waves  propagating  from  a homogeneous  region  into 
a nonuniform-depth  regime.  The  spectral  transformation  at  a point 
X along  a ray  that  originates  at  in  the  homogeneous  region  where 
ii  “ is  described  via  (42)  and  (46  a,b),  e.g.,  on  setting 

s=0  where  x = x for  convenience,  we  have 

— - - P3  ' 


ij<^°\k,x)  = ) 


x(s)  = X + / k(x. 

— —CO  / — — 

-'0 


k(x)  ds 


k(x)  = k 


/C  } ds 
g 


(49  a,b,c) 
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The  solution  to  the  spatial  transformation  of  (_k,2c)  via  ray 

training  can  involve  excessive  numerical  computation  In  particular  in 
a situation  wlilch  requires  the  determination  of  for  several  points 

in  the  inhomogeneous  region.  An  alternate  approach  is  to  solve  (A5) 
together  with  (8  a)  (l.e,,  ^x_k  = 0),  using  a proper  finite  difference 
scheme  in  the  x-plane.  The  third  and  possibly  best  alternative  is  to 
take  advantage  of  both  of  these  and  use  the  fact  that  the  transformation 
of  really  involves  only  the  spatial  modulation  of  Ic-space.  Therefore, 

the  whole  problem  of  predicting  (k,x)  reduces  to  tlie  prediction  of  a 

Ic-space  in  the  inhomogeneous  region  corresponding  to  a specified  in  the 

homogeneous  region,  l.e.,  to  the  solution  of  ^x_k  = 0 subject  to  the 
condition  Jk  = along  the  homogeneous  boundary.  This  tlien  enables  us 
to  set  ii)^^^(k,x)  = simply  because  that  particular  solution  to 

^Xjc  = 0 implies  that  we  can  trace  a wave  ray  starting  with  Jc  at  point  jc 
in  the  inhomogeneous  region  back  into  the  homogeneous  boundary  where 
_k  1^.  This  approach  is  simplified  further  by  working  in  a co, 0-space  rather 
than  In  _k-space  since,  in  this  case,  lo  is  spatially  invariant  by  virtue  of 
(8  b)  and,  therefore,  it  would  suffice  to  predict  only  the  spatial  modula- 
tion of  0 In  terms  of  a specified  homogeneous  boundary  value  0^.  In  doing 
so,  it  is  advantageous  to  introduce  the  spectral  density  n^^^(m,0,2c)  over 
a oa, 0-plane  defined  as 

= k 1^  (50) 

so  that  dud0  = Now,  wo  can  easily  show,  using  Ic  = k(cos  0,  sin  0) 

in  ^xjc  “ 0,  expanding  and  utilizing  the  dispersion  relation  (17)  for  (3k/3x) 


23 


and  (3k/8y)  in  the  resulting  expansion,  that 


c • ve 

_g  _ 


U)  " 

slnh  2q  — 


Vh 


(51) 


with  n = (-  sin  9,  cos  9),  Is  the  required  relation  governing  the  spatial 

modulation  of  9.  Therefore,  (51)  with  9=9  being  prescribed  along  the 

00 

homogeneous  boundary  and  the  equation  resulting  from  substitutions  of  (50) 
and  (33)  into  (42),  i.e.. 


s/°^(a),9,x)  = ^ 


nf°\(o,9  j 


C52) 


are  sufficient  for  the  spectral  prediction  problem.  The  following  example 
will  serve  to  illustrate  this  prediction  for  random  waves  propagating  from 
deep  water  into  shallow  water. 


SPECTRAL  PREDICTION  IN  SHALLOW  WATER 


Consider  as  an  example  an  underwater  topography  illustrated  in  Figure  2. 
We  assume  that  waves  propagate  from  deep  water  (i.e.,  x = -“)  over  an 
intermediate  region  with  parallel  isobaths  into  shallow  water  where  the 
variation  of  the  still  water  depth  is  described  (in  meters)  by 


h(x,y) 


2{1  - 10“^  x^Il 


g-l/2(y/10)2^j 


X > 0 , |y|  < “ (53) 


X (METERS) 


2 
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The  homogeneous  deep  water  wave  field  is  specified  in  terms  of  a spectral 
density  of  the  form 


where 


(U),0  J 

CO 


f (9  J f (“) 

1 2 


(54) 


f (0  ) _ 8 

I 00  =S  ■ ~ 

1 3Tr 


4 

cos 


0 

00 


|0  I < tt/2 
' 00  ' 


(55)  ' 


represents  a directional  spreading  factor,  and 

= (0.0081)  g^oj"^  exp[-  0.74(|^)'^J  , u > 0 (56) 

is  the  well-known  Pierson-Moskowitz  spectral  form,  with  U representing  the 
overwater  wind  speed.  We  take  U = 1.5g  as  a specific  case  here.  The 
corresponding  directional  spectral  density,  (w»6^)  > is  illustrated  in 

Figure  3 showing  various  contours  of  = constant. 

In  the  region  from  deep  water  to  x = 0 where  h = 2m,  the  wave  field 
is  inhomogeneous  only  in  the  x-direction.  Therefore,  (51)  or,  more 
simply,  ^xk  = 0,  Implies  that 

k sin  0 = k sin  0 = const  , (57) 

00  00 

which  is  Snell's  law  of  refraction  appropriate  to  an  underwater  topography 
with  parallel  Isobaths.  Therefore,  it  follows  from  (52)  and  (57)  with 
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CC  ) = u)/2k  that 

g OO  tiO 


>(0), 


k.  (jJ 


il  C“»0i2i)  “ 2 'k 

“ g 


sin  sin  0|1  , x < 0,  | sin  0|  < 1 


(58) 


For  X ^ 0,  the  isobaths  are  no  longer  parallel.  In  this  region,  the 

significant  spectral  components  of  the  wave  field,  e.g.,  those  for  which 

2 

^ — 1.2,  are  in  shallow  water,  i.e.,  (u  h/g)  < (ir/lO).  Further,  it 
is  evident  Irom  (57)  that  propagation  directions  of  these  shallow  water 
components  are  very  much  focused  towards  the  x-axis.  These  conditions  will 
enable  us  to  derive  an  approximate  closed  form  solution  for  x > 0 as  follows. 
First,  we  expand  (57)  in  the  form 


90. „ .«90  w 

C cos  0 — + C sin  0 — = . „ 

g 9x  g 9y  sinh  2q 


(sin  0 - cos  0 ■^) 

9x  9y 


(39) 


Using  small  angle  approximations,  i.e.,  sin  0=0  and  cos  0=1,  and  keeping 
terms  of  the  first  order  in  0,  (59)  becomes 


C 


90 


g 9x  sinh  2q  9x 


From  the  definitions  of  phase  speed 


1^9  = 


9h 


sinh  2q  9y 


(60) 


C = w/k  = (q  tanh  q)^^^/k 


and  Cg  given  by  (33),  it  follows  that  (60)  simplifies  to 


(61) 


28 


= i’i'. 

Sx  C 3x  C 3y 


(82) 


or,  simply 


c-^}  = 


-c-2^ 

<iy 


(63) 


Integrating,  we  will  obtain 


0(x,y)  = C(x 


,y)  [-/ 


C-2  3£  Jx  + » I 
Q 3y  C(0,y)  f 


(64) 


Since  C = [g  h(x,y)]^^^  in  shallow  water,  we  can  rev;rite  (64)  as 


a/  ■>  1 ( f -3/2  Oh  .w  + 2_l(P^y)_  I 

8<x,y)  . 2 I,  (K,y)  I.  - ..=■  i 


(6  .) 


For  the  example  of  interest  here,  with  h cleaned  by  (53),  this  becomes 


0(x,y)  = byd  + e 


by2j-l  ( 


h(x,y) 


a(l  + e ) 


, -Ir,  h(x,y),l/2  ) 

sin  [1 - X 


-3  -2 
where  a = 2 x 10  and  b = 0.5  x 10 

Figure  4 illustrates  the  spectral  transformations  at  three  points 
l^(14.14ra,  0)  , P2(15.78m,  10m)  and  F^(18.76ra,  20mO  of  Figure  2,  where 
h ■=•  1.2m.  It  can  be  shown  from  (57)  and  (66)  tiiat 
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6(14,14,0)  = 0.775  0(0, y)  = 0.775  sin"  {—  ;nii  0 I 

k 


0(15.76,  10)  = -0.048  -f-  0(14.14,0) 


0(18.76,20)  = -0.036  + 0(]4.14,O) 


1^,  P2  and  F^,  respectively.  Consequently,  the  d.-iisitic-s  ;it  the:'.:  iluei. 


points  are  identical  except  for  a uniform  shift  of  approxiit-i t ■ 


I ■ t T 


the  directions  of  spectral  components  at  1=  , and  a shift  of  aiqe  ...r  i - 
-2.06“  for  those  at  relative  to  the  corresponding  components  at  P.  . 


REAL  REPRESENTATIONS 


In  order  to  utilize  results  such  as  those  given  by  (25  n,b)  ree.di  .y, 
it  is  appropriate  to  express  them  in  terras  of  real-valued  quant! tie".  Foi 
a real-valued  random  process  qCl^jt), 


{dA^^^  (k,x)  }*  = dA^*^^(-k,x) 


It  follows  then  that  we  can  define  two  real  random  processes  and 

such  that 


dU^°^  = dA*^°^  + {dA^°h* 


(68  a,h) 


It  is  evident  that  is  even  whereas  is  odd,  i.e.. 
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dU^°\-k.,x)  =■  dU^^^Ck.x) 


dV^°\-k,x)  = -dV^^^Qc  x) 


Furthermore,  we  can  show 


(69  a.b) 


E{dU^^^  dV^^h  = 0 


E{[dU^°^]2}=  E{[dV<^°^]2}=  2 E{|dA^°^|  } = 2 i);^°\k,x)  6k 


Using  C68  a,b),  we  can  rewrite  dA^  ' in  the  form 


dA^°^  = i {dU^°^  - 1 dV''°h 


Making  use  of  (71)  and  also  the  fact  that 


F(-k,x,z)  = - F(k,x,z) 


one  can  transform  all  the  results  previously  given  in  a complex  form  into 
real  representations.  In  particular,  one  can  verify  that  the  free  surface 
(25  a)  becomes  to  the  order  0(|Vh|^) 


n(x.t)  = ^ J [ 


F dV^®^}  cos  X + + F dU^^^}  sin  x] 


This  can  be  further  simplified  if  we  express  dA'  ' in  terms  of  its  amplitude 
jdA^^^I  and  phase  y,  i.e.. 
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1 


(lA 


(0) 


= <1A 


(0) 


Correspondingly,  one  would  get 


. 2 Ua<“>|  CO. 


= -2  jdA^^^l  oil 


On  substitution  from  (75  a,b),  (73)  becomes 

ri(x,  t)  = {cosCx  + ^')  ^ 

k 

It  should  be  emphasized  that  F In  the.  precedin  . -.  .nr'  :--:!  ■ . , 

F “ F(_k,x,  0)  = q tanh  ■;  c,  . 

where  a^,  a^,  are  as  provloiusly  detined  by  (21  a,:-  ■ ' ' ■■  ‘ ■■  by 

. This  function  incorporates  nonuniforin  deptt.  f.l  ■ 

OCl^h!),  and  constitutes  a modification  to  the  r . -.a  .'i;  i n 

represented  by  the  cosine  term  in  (76).  Hoxv-ever,  ■•’.s  pr*'\'-  asi  % 

the  first  two  order  solutions  are  orthogonal  so  tn.'-t^  al; 

surface  is  modulated  by  the  nonuniform  depth  t..  ■ li  i d!  r x ' • > 

no  modification  to  the  free  surface  spectral  di  o i 

To  examine  the  effect  of  a sloping  bottom  ' .■  i . 

explicitly,  consider  the  case  of  parallel  bottui.i  • ■ .mi  ,:i  ' ; 0, 

3h/9x  ^ 0.  Under  these  '.onditlons,  it;  can  be  vcrift'.  fr-  "c), 

(57)  and  (A. 4 a,b,c)  that  (77)  simplifies  to 


J 
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F(k,x,0) 


q ' (3  slnh  2q  + 2q)  dh 


(slnh  2q  + 2q)' 


dx 


cos  0 


(78) 


This  is  illustrated  in  Figure  5,  which  shows  that  the  effect  of  a sloping 
bottom  initially  increases  at  first  as  waves  feel  the  bottom,  then 
monotonically  decreases  at  shallower  depths,  both  effects  being  most 
pronounced  for  components  at  normal  incidence  (i.e.,  0 = 0°)  to  the  bottom 
contours  as  expected.  Consider  now  the  net  effect  of  this  modulation  to 
the  locally  flat  bottom  solution  for  the  free  surface  regarded  as  a sum 
of  constituent  wavelets  of  the  form 


cos(x  + p)  + F sia(x  + p) 


(79) 


It  is  evident  that  F has  the  same  sign  as  the  local  dh/dx.  Hence,  when  dh/dx  < 0, 
as  we  move  from  a crest  Ce-8«»  x ■*"  P “ 0)  of  an  elementary  component  cos(x  + p) 
to  a following  trough  (x  + U “ if)#  F sinCx  p)  <0.  On  the  other  hand,  as  we 
move  up  from  the  same  trough  up  to  the  next  crest,  i.e.,  it  < x + P < 2tt,  we 
see  that  F sin(x  + p)  >0.  It  follows,  therefore,  that  the  profile  of  an 
elementary  component,  cos(x  + p),  of  the  locally  flat  bottom  solution  becomes 
assymmetrical  by  a suppression  from  a crest  to  the  trough  followed  by  a magnifi- 
cation from  the  trough  to  the  next  crest.  Clearly,  the  same  profile  remains 
unaffected  at  any  crest  or  trough  since  sin  (x  + p)  = 0 at  such  points.  The 
situation  is  exactly  reversed  when  dh/dx  > 0. 

PROBABILITY  STRUCTURE  AND  ENVELOPE  PROCESS 


The  probability  structure  of  a wave  field  sucli  as  that  depicted  by  (76) 
is  embedded  in  the  random  character  of  the  complex  differential  amplitude 


FICTRli  5.  Ami) ).i tilde  Modulation  of  First  Ordet 
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in  general,  in  terms  of  both  its  amplitude  jdA^^^I  and  phase  p. 
However,  any  rational  specification  as  to  the  explicit  forms  of  these  must 
be  guided  by  the  probability  structure  of  Ti(x»t)  itself.  In  this  case  it 
is  widely  accepted  that  this  structure  is  Gaussian.  This,  therefore,  leads 
to  a model  of  random  phases  for  the  amplitude  process  dA^*^^  which  is 
well-known  in  theory  and  applications  in  the  homogeneous  cases.  The  model 
of  random  phases  simply  means  -as  the  name  implies-  that  |dA^^^|  is  a 
deterministic  function  whereas  the  phases  p are  uniformly  distributed  over 
CO,  2ti)  and,  therefore.  Independent.  In  the  homogeneous  case  both  |dA^^^| 
and  p are  simply  functions  of  Jc  or,  equivalently  u and  9,  without  any  depen- 
dence on  X.  For  example,  in  deep  water 


dA^^^(k  ) = ) 6k 


provides  the  explicit  form  required  so  that 

E{ldA^°^|^}  = ) 6k  C81) 

An  obvious  representation  for  dA^*^^(jc,x)  in  the  spatially  inhomogeneous  case 
is  in  the  same  form,  i.e.. 


dA^°^(k.2i)  = [i(^^°\k,x)6k]^^^ 


where,  for  each  Jt,  p(]t,x)  is  distributed  uniformly  over  (0,  2tt)  . For  a 
wave  field  which  originates  from  a homogeneous  region,  (35)  implies  that 


r 


M(k,x)  = P(]S„) 


as  a general  Invariant  property.  In  this  manner  one  can  obtain  from  (76),  on 
substitution  from  (82),  a spatially  Inhomogeneous  random  wave  field  which 
is  Gaussian  in  probability  structure  [see.  Appendix  C],  l.e., 


f^(u;x)  clu  = Prob{u  j<  ^(ic,  t)  ^ u + du} 


1 r 1 , u. 2,  j 

exp  1-  ^ 

n 


where 


o ^ = o ^(x)  = E {ri^(x,t)) 
n T)  — — 


One  of  the  v\seful  descriptions  associated  with  a random  wave  surface 

of  the  form  (25  a)  is  its  envelope  process  which  is  defined,  following 

Cramer  and  Leadbetter  [1967],  at  a fixed  x as 

~o 


:(t)  = [n“(x  ,t)  + n^C^:  ,t)]' 

— O — O 


where  S(t)  > lri(x  ,t)j  and  iE(t)  = |n(x  ,t)|  when  n(^  ,t)  = 0.  It  is 
known  [see,  e.g,,  Crandall,  1970;  Yang,  1972]  that  tl\e  envelope  definition 
of  (86)  is  essentially  the  same  as  that  due  to  Rice  [1955]  for  narrow- 
band  processes.  In  (86),  r)(x^,t),  whicli  is  called  the  Hilbert  cransfona  of 

A 

n(2^,t),  is  obtained  by  passing  through  a filter  with  the  frequency 

response 


P 
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G((ij) 


1 , U)  < 0 


0 , 0)  = 0 


-i  , u > 0 


(87) 


or,  simply  G(u))  = -i  sgnCw),  where  sgnCw)  designates  the  sign  of  w.  Therefore, 
it  follows  i-om  (87)  and  (25  a)  that 


h(>EQ.t) = -i  y 


sgn(w)  [1  - iF(k.x  dA^^^ (k,x  ) exp(ix)  (88) 


Now,  by  making  use  of  (67)  through  (72),  and  associating  -k  with  -m  and  +k 
with  +0),  we  can  rewrite  (88)  in  the  following  equivalent  real  form: 


n(x  ,t)  = J'  sgniui)  |dA^^^  I {sin(x  + y)  + F cos(x  + y)) 


(89) 


Assume  now  that  n(x  ,t)  is  Gaussian.  It  follows  then  that  n(x,t)  is 
— o 

also  Gaussian.  Furthermore,  we  can  verify,  from  (76)  and  (89)  with  y being 
uniformly  distributed  in  (0,  2it)  , that 


o~^(x  ) = a ^(x  ) = E{n^(x  .t)} 
r|  — o n ~~o  —o 


E{n(2SQ,t)  n(x^,t  + t)}  = E{ri(j^,t)  n(x^,t  + t)} 

= 4 J' ) (cos  toT  + F sin  wt}  6k 
k 

E{n(x^,t)  Ti(5^,t)}  = ysgn((o)  F(k,2^,0)  (k,2^^) 


6k 


(90  a,b,c) 
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Letting, 


P(x„)  = 


E{ri(x^,t)  nCx^.t)} 


the  joint  probability  density  of  n(x  ,t)  and  ri(x  ,t)  becomes 

— o 


2.1/2 

2v  a (1  - p ) 

n 


exp  {- 


2 « , ''2 

n - 2p  nn  + n 
^ 

0 2,,  2,  ' 

2 o (1  - p ) 

n 


< n,n  < 


Now,  we  can  use  the  fact 


0 ''2  T / ? 

Prob{  H(t)  < u)  = Prob{[n‘'  + n ] < u} 


to  show  that  the  density  of  HCt)  is 


, / ■,  u 1 u I 

- ^ ^ ; — 2'n  " 7^;, — X I 

2ti  0^  (i  - P ) I 2 (1  - p ) ; 


i i’  sin  2v  I , 

J I 2 T i 

I 2 o ^1  - p^)  ) 

n 


Recognizing  further  that 
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where  I [•]  represents  the  zero-order  modified  Bessel  function  of  the  first 
o 

kind,  (S3)  becomes 


Cu;Xo) 


a - p ) 


2 o ^(1  - p^)  I 

h 


2 2 ’ 
2 0 ^(1  - p^)  J 


0 ^ u < “ (95) 

Finally,  it  can  be  verified  that  the  j-th  moment  of  the  envelope  process 
is  in  the  form 

E{  = TT"^^^(2a^)^Cl  - p^)  ^ r(i  + 1)  r(j  + |)  F(^  + 4 + r 

(96) 

where  F and  ^ are  the  gamma  function  and  the  hypergoemetric  function,  respec- 
tively. The  distribution  of  wave  heights,  H = 2H  , follows  immediately  from 
(95),  i.e., 

f ■ i f , («) 

2 ^2  2 

Hence,  by  recognizing  that  E{H  } = H = 4 , we  obtain 


i„(v;x  ) = Y J _ 

ir(l  - p*-)  2H  (L  - p -)'  21i^(l  - p^). 


0 < V 


The  only  constraint  on  the  preceding  results  Is  tlie  assumption  that 
the  free  surface  ri(x,t)  Is  Oaussian.  In  deep  water  or  in  a region  of 
uniform  depth  where  F 0,  p = 0 and,  therefore,  (98)  and  (99)  become 
the  Rayleigh  density  and  the  associated  j-tli  moment,  respectively,  i.e., 


V 1 V 

exp< 


2ld’  ’ 


0 < V < 


ECH^]  = (2  11'^)^''^  r(l; 


i-  1) 


(100  a,b) 


These  results,  which  are  widely  used  in  theory  and  applications  with  the 
inherent  constraint  that  the  underlying  free  surface  spectrum  is  narrow- 
band,  are  in  fact  valid  without  any  bandwidth,  restrictions,  it  is,  there- 
fore, not  surprising  that  most  empirical  wave  height  frequency  curves 
obtained  in  field  applications  involving  wide-hand  spectra  support  this 
conclusion.  Nonuniform  depth  effects  on  the  distribution  of  wave  heights 
are  characterized  to  the  order  0([vh|“)  l;y  p,  which  is  simply  the  correlation 


coefficient  of  n(x  ,t)  and  n(x  ,t).  Tills  correlation  Introduces  a skewness 
— o — o 

in  the  Rayleigh  form  appropriate  to  a flat  bottom  towards  lower  wave  height 
values.  For  example,  It  can  be  verified  frc’m  (99)  and  (100  b)  with  J “ 1 that 
the  ratio  between  the  mean  values  of  H corresponding  to  a variable  bottom 


2 2 

and  a flat  one  is  (1  - P )/^(5/4;  3/4;  1;  p ),  whlcli  is  always  less  tlian 
2 

unity  since  p <1. 


CONCLUDING  REMAliKS 


Propagation  and  refraction  of  linear  waves  in  water  of  nonuniform 
depth  have  been  extended  to  the  case  of  random  waves  characterized  by  a 
two-dimensional  spectral  distribution.  Attention  was  focused  to  the 
orthogonal  representation  of  the  free  surface,  its  spectral  density  and  to 
the  statistical  distribution  of  wave  heights  with  the  primary  objective 
to  incorporate  nonuniform  depth  effects  explicitly  into  well-known 
solutions  associated  with  a locally  flat  bottom.  Although  the  methodology 
employed  corresponds  exactly  to  the  linear  approximation  of  that  developed 
by  Chu  and  Mei  [1970]  for  the  case  of  monochromatic  waves,  the  solution 
derived  here  for  the  order  0(|^h|)  is  different,  and  believed  to  be  the 
valid  one.  What  is,  however,  more  important  is  the  extension  of  the  results 
to  a random  wave  field  and  to  a continuous  spectral  representation.  This, 
therefore,  compliments  and  generalizes  the  deterministic  concepts,  leading 
to  various  probabilistic  descriptions  of  wave  properties  in  spatially 
inhomogeneous  media.  The  probabilistic  description  was  also  limited  to  the 
free  surface  and  its  envelope.  The  extension  of  these  to  the  surface  gradient 
and  other  kinematic  and  dynamic  properties  such  as  pressure,  velocities,  etc. 
does  not  offer  any  special  difficulty  and  needs  to  be  done.  In  particular, 
the  predictions  of  wave  breaking  and  breaker  line  for  waves  Impinging  on 
a beach  can  be  formulated  as  a smooth  generalization  to  what  Biesel  [1952] 
and  Gaughan  and  Komar  [1973]  have  done  in  the  case  of  one-dimensional 


monochromatic  waves.  It  must,  however,  be  pointed  out  that  a definite 


r 
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tiiscr epancy  exists  between  tiie  solntion  of  the  order  O(l^lij)  here  and  the 
corresponding  one  used  by  these  authors.  These  and  various  other  aspects 
of  the  problem  such  as  an  extension  to  Include  nonli neari ties  ai e considered 
for  future  study. 
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APPENDIX  A - Derivation  of  Coefficient  dC 


(1) 


To  determine  we  use  the  condition  that  as  q dB^^^  = 0; 


Therefore,  from  (19)  and  (18  b) , 


dC 


(1)  _ )(g/“)  dA^"^  F(k,i,E)) 


(0) 


lim 


(A.l) 


cosh  q 


or,  on  substitution  from  (20)  and  (21  a,b,c). 


(u)/g)  dC^^^  = lira  dA^°^{a  ^ — + ~ ^ ^ + a.  } (A. 2) 

^ cosh  q 2 cosh  q 3 cosh  q 

q->-  oo 


We  can  now  verify,  by  using  (lA  a,b)  and  (18  b)  in  (16)  and  integrating. 


that 


k 

k 


V dA 


(0) 


dA 


(0) 


k 2 V q (q  tanh  q - 1)  ^k 
k sinh  2q  + 2q  2k^ 


(A.3) 


where,  from  (17), 


Vk 


2k^  jh 

sinh  2q  + 2q 


V 


k 


2k 


3 


sinh  2q  + 2q 


Ih  + i 
^^2 


Vq 


= k ^h(- 


sinh  2q 
sinh  2q  + 2q 


■) 


(A. A a,b,c) 


f 


h5 


On  the  basis  of  the  preceding  and  using  the  definitions  of  5^,  and 

5^  given  by  C21  a,b,c)  it  immediately  follows  that  dC^^^  = 0. 
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APPENDIX  B - Coordinate  Transformation  and  Derivation  of  Equation  (32) 

Consider  a change  of  variables  from  x = (x,y)  to  an  orthogonal  curvilinear 
system  Cs,n),  where  s and  n denote  coordinates  measured  along  and  perpendicular 
to  the  orthogonal  trajectories  (rays)  of  a component  wavelet,  as  schematically 
illustrated  in  Figure  1.  Referring  to  the  same  figure,  we  see  that 

- ^ 

_k  = = Ceos  0,  sin  6) 

= (-  sin  0,  cos  0)  (B.l  a,b) 


represent  unit  vectors  in  the  direction  of  s and  n at  x.  Furthermore,  desig- 
nating the  separation  distance  between  two  neighboring  s trajectories  by 
p,  we  have 


ii  = i iE 
3n  p 3s 


Now,  for  an  arbitrary  function  f = f(x,y),  we  can  write 


9 s 


k 


Vf 


3 f 

cos  0 -z—  + sin  0 

dX 


ll 

3y 


(B.2) 


= n 
3n  - 


Vf 


sin  0 + cos  0 

3x 


3y 


(B.3  a,b)  j 

I 

1 


47 


, on  solving  for  3f/8x  and  3f/3y 


3t 

3x 


cos 


sin 


0 


3n 


3y 


sin 


cos 


We  can  then  use  (B.4  a,b)  for  cos  0 and  sin  0,  respectively, 
obtain 


(cos  0)  = cos  0 (cos  0)  - sin  0 ~ (cos  0) 

(sin  0)  = sin  0 — (sin  0)  + cos  0 (sin  0) 


Combining  (B.5  a)  and  (B.5  b) , and  utilizing  (B.2)  yields 


V 


k = V 


i l£ 

p 3s 


The  group  velocity  is  defined  by 


C 

-g 


= C k = 
g - 


0)  /Sinh  2q  + 2q. 
2k  sinh  2q  ' 


k 


(B.4  a,b) 


to 


(.B.5  a,b) 


(B.6) 


(B.7) 


In  deep  water  as  q this  becomes  a constant  Independent  of  x,  i.e., 


(B.8) 


Now,  using  the  definitions  (B.7)  and  (B.8),  we  can  show  that 


— — fin  = 2(q  tanh  q - 1) 

2 8q  C sinh  2q  + 2q 

O 


(B.9) 


On  the  basis  of  the  preceding  results  and  definitions,  (31)  becomes 


dA 


(0)  3s  ^ 


= la  - -1 

2 3q  Cg  ^ 3s  2p  3s 


(B.IO) 


or,  simply 


^ _i  {dA^°)}  = 1 -i.  {in^^S^} t lE 


dA 


(0)  3s 


2 3s 


2p  3s 


(B.ll) 


Realizing  that  (C  )^  - constant  for  a fixed  k,  and  Integrating  (B.ll), 

O 

we  obtain  the  desired  result  that 


^1/2  „ 1/2  .,(0),,  , 

p C dA  (k,x)  = constant 

O 


(B.12) 


along  a wave  ray. 


APPENDIX  C - Proof  of  Gaussian  Proper! 


A mean-zero  random  process  is  Gaussian  if  its  log-characteristic  function 
possesses  the  following  form: 

In  - - I // / 2-2-  'r  '^2'  - 


X X T T 


t(x2»  ^2)  dXj^  dx2  dt^  dt2 


(C.l) 


where 


^nn^-i’  -2’  *^i’  " E{nCX]^,  n*(x2»  <^2)} 


is  the  covariance  function  of  ri(Xft).  Recall  from  (76)  that  we  have,  to 


0(1  Vhp) 


n(x,t)  = / 


|dA  I {cos  (x  + p)  + F sin  (x  + P) ) 


CC.3) 


Hence,  it  follows  from  (C.3),  (C.2)  and  the  fact  that  p*s  are  independent 
random  variables  uniformly  distributed  in  (0,  2::)  that 


-2’  *^l’  *^2^  " *^nn^-l’  -2’  ^2  ~ ‘^1^ 


(C.4) 


■|  /|dA{°^|  |dA2^°^|  (cos  (x^  - X2)  + - F2)  sin  (x^  " X2)  > 

k 


where 


IcIaJ^^I  = |dA^°\k,x^)| 


= IdA^^^k.xpl 


f-2 

- X2  = - V k • d X - toCt^  - t^)) 


Fi  - F2  = F(k,x^,0)  - F(k,X2.0) 


(C.5  a,b,c,d) 


By  definition 


ff 

J J nOf.  I 


M {T(x,t)}  = Efexp  [i  J J nOf.t)  T(x,t)  dx  dt} 

X T 


(C.6) 


Now,  we  approximate  (C.3)  by  a Riemann  sum: 


flCx.t)  = I |dA^^^|{cos  (x  +p  ) + F sin  (x  +P  )} 
jj  jji  nm  ' nm  nm  nm  ^nm  nm 


where 
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F = F(k  , X,  0) 
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= p(k  , x) 
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(C.8  a,b,c,d) 
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From  CC.6),  (C.7)  and  the  fact  that  M are  uniformly  distributed  in 

nm  ■' 

C0,2'ir),  it  is  immediate  that 


M (rCx.t)}  = II  z ^ 
n ““  nm 

n,m 


(C.9) 


where 


Zir 

„ , exp{i(a  cos  li  - b sin  p )}  dp 

nm  Ztt  / Q nm  nm  nm  nm  nm 


■ “ fa 
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X T 
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^ ~ J J {cos  X ™ + F sin  \ } xCx.t)  dx  dt 

nm  tr  m nm  ' ''nm  nm  ^nm  — ' — 


,(0) 


J J I'^^nm  - ^nm  ^nm^ 

X X 


Equation  (C.IO  a)  can  be  integrated  as  follows: 


2tt 

Z = “ / exp{i  Y cos(p  - 6 ) } dp  = J (y  ) 

nm  2tt  / „ nm  nm  nm  nm  o nm 


(C.ll) 


where  is  the  zero-order  Bessel  function  and 


Y = [a^  + b^ 

nm  nm  nm 


- ^ “1  / nm. 

6 = tan  C-  ) 

nm  a 

nm 


(C.12  a) 
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Therefore,  it  follows  CC.9)  and  (C.ll)  that 


In  M {T(x,t))  = I In  J (y  ) 
n — o nm 

n,m 


(C.13) 


Now,  If  we  let  = max  Ak  , which  can  be  made  as  small  as  we  please 

— nm  —run 

X 


by  choosing  sufficiently  small,  it  becomes  evident  from  (C.8  a)  and 


1/2. 


(C.IO  b,c)  that  a and  b would  be  G(Ak  ) at  the  most.  Therefore, 
^ ' nm  nm  — nm 


,l/2s 


Y , which  is  of  the  same  order  0(Ak  )by  virtue  of  (C,12  a),  can  be 
nm  ~nm 


made  sufficiently  small,  enabling  one  to  vrrite 


J(y  )=l--fY^  + 0(y^  ) 

o nm  4 nm  nm 


(C.IA) 


In  this  manner,  we  can  rewrite  CC.13)  as 


In  M {T(x,t)}  = I ln{l  - j + 0Cy^*,.)J 
n ^ 4 ' nm  nm 
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= I {-f  Y^  + OCy"*  )} 

^ 4 nm  run 
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(C.15) 


where  the  expansion 


ln{l  - X Y^  } 
4 nm 


- T + 0(y^^) 
4 ' nm  nm 


(C.16) 


has  been  used. 


With  the  substitution  of  Y » as  defined  through  (C.12  a),  and  (C.IO  b,c), 

nm 


into  (C.15),  the  log-characteristic  function  can  be  written  as 
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In  M {t} 
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(C.17) 

Note  that  as  max  “ max  /Me  -♦■  0,  the  square-bracketed  double  sum  in 

n,m  Xj*^*”* 

CC.17)  approaches  2 K^^(2Cj^,  3C2»  tj^,  t2)  and  (C.17)  becomes  Identical  with 
(C.l).  Therefore,  n(x, t)  is  a Gaussian  random  process. 
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RANDOM  WAVE-CURRENT  INTERACTIONS  IN 
WATER  OF  VARYING  DEPTH 

M.  A.  Tayfun,  R.  a.  Dai.rymple  and  C.  Y.  Yang* 

IJeparinient  of  C ivil  Engineering,  University  of  lielaware,  Newark,  Delawfarc  19711,  U.S.A. 

.Abstract  Refraction  of  incoherent  random  gravity  waves  with  currents  and  bottom  topo- 
graphy results  in  spatial  variations  in  the  spectral  characteristics  of  the  free  surface.  Prediction 
of  such  variations  ha.scd  on  the  radiation  transfer  equation  is  in  a simple  analytic  form  for  the 
case  of  one  dimensional  inhomogeneities  in  currents  and  topography.  This  analytic  form  is 
esamined  in  terms  of  two-dimensional  wave  number-  and  polar  frequenej-  direction  spectra 
along  the  associated  dynamic  and  kinematic  constraints  relevant  to  w'ave  breaking  and 
retlection.  Results  are  specialised  to  the  simplest  case  of  horizontal  shear  currents  in  deep  and 
shallow  water  w ith  explicit  examples  to  illustrate  the  relative  and  combined  elTects  of  currents 
and  topography  on  free  surface  spectra. 

INTRODUCTION 

Refraction  of  surface  waves  interacting  with  currents  and/or  underwater  topography 
results  in  spatial  variations  in  their  kinematic  and  dynamic  properties.  At  present,  a well- 
known  wave  refraction  theory  exists  to  predict  such  variations  for  linear  monochromatic 
waves.  However,  description  of  the  sea  surface  in  terms  of  monochromatic  waves  is  at  best 
an  approximation  since  most  realistic  sea  states  have  a more  complex,  randomly  irregular 
structure.  This  structure  is  concisely  characterized  by  a two-dimensional  mean  square 
spectral  distribution  over  a wave  number  space  or  a polar  frequency  direction  space. 
Consequently,  the  application  of  the  concepts  developed  for  moniKhromatic  waves  inter- 
acting w ith  currents  and  underwater  topography  to  the  ptediction  of  the  spatial  transforma- 
tion of  spectral  characteristics  is  of  interest  in  the  study  of  random  waves  and  related 
phenomena.  The  earliest  effort  in  this  direction  is  by  Longuet-Higgins  (1956,  1957),  who 
derived  the  transformation  of  two-dimensional  wave  number  spectra  by  refraction  over  a 
general  topography.  Later,  Phillips  (1966),  Hasselmann  (1968)  and  others  generalized 
Longuet-Higgins’  results  as  a systematic  energy  balance  formulation  for  the  piediction  of 
spatially  inhomogeneous  spectra,  taking  into  account  current  interactions,  various  dissipa- 
tive and  generative  effects  as  well.  This  formulation  has  been  recently  extended  to  nonlinear 
random  waves  by  Willebrand  (1975).  However,  explicit  solutions  and  applications  of  the 
energy  balance  or  radiation  transfer  equation  have  been  demonstrated  only  in  a few  cases. 
Karlsson’s  (1969)  numerical  predictions  of  the  refractive  transformation  of  polar  frequency- 
direction  spectra  over  topographies  with  parallel  and  irregular  contours,  Krasitskiy's 
(1974)  closed  form  predictions  of  a similar  nature  but  restricted  to  parallel  contours  in 
essence  constitute  further  explorations  and  applications  on  Longuet-Higgins’  (1957)  depth- 
refraction  solution.  Similarly,  by  also  taking  into  account  an  approximate  form  of  bottom 
friction,  Collins  (1972)  numerically  computed  spatial  variations  in  spectra  due  to  depth 
refraction.  In  the  case  of  unidirectional  waves  interacting  with  an  opposing  or  following 

•Please  direct  ail  correspondence  to  C.  Y.  Yang. 

40.^ 


r 1 


■104  M.  A.  R a,  Dairvmi'ii  and  t . Y.  Yanc. 

current  in  deep  water,  the  explicit  translorniation  of  one-dimensional  frequency  spectra 
was  given  by  Phillips  (1966,  p.  60).  The  applications  of  this  to  current  measurements  and 
wave  forces  were  demonstrated  by  Huang  ci  al.  (1972)  and  Tung  and  Huang  (1973). 

The  motivation  here  is  to  explore  some  straightforward  solutions  of  the  radiation  transfer 
equation  for  a medium  with  one-dimensional  variations  in  horizontal  current  components 
and  topography.  Attention  is  also  restricted  to  linear  waves,  neglecting  various  nonlinear 
mechanisms,  dissipation  and  generation  elTects.  Guided  by  analogous  concepts  on  mono- 
chromatic waves,  dynamic  and  kinematic  constraints  related  to  wave  breakingandreflectionts 
are  examined.  Results  are  specialized  in  detail  to  the  simplest  case  of  horizontal  shear  currents 
in  deep  and  shallow  water,  and  their  salient  features  are  illustrated  with  explicit  examples. 


DEFINITIONS  AND  ANALYSIS 
An  incoherent  random  wave  field  can  be  represented  by 

q (x,  I)  ---  I a„  (x,  /)  cos  (k„  . X — on  |x„)  (1) 

where  ix„  represents  random  phases  uniformly  distributed  over  (0,  2n);  k„  and  io„  define 
the  vector  wave  number  and  frequency,  respectively;  and  a„  is  a Fourier  amplitude  regarded 
as  a slowly  varying  function  of  time,  I,  and  position,  x = (.v,  >)  in  a horizontal  co-ordinate 
system  fixed  at  still  water  level.  As  the  vector  wave  numbers,  k„,  are  densely  distributed 
over  a k (A’l,  k.^)  plane,  the  corresponding  amplitudes,  a„,  tend  to  become  dense  and  also 
infinitesimal.  Under  this  condition,  the  function,  y,  represented  by 

V (k;  X,  f)  dk  ^ ^ o;,  (2) 

with  dk  dA,  dA..,  as  a shorthand  notation,  is  defined  as  the  wave  number  spectral  density, 
describing  the  distribution  of  mean  square  surface  deformation  over  k-space  locally,  i.e. 

■ q- (x,  t)  ^ o;,  I v(k;x,  Odk.  (3) 

2 " J 


I his  distribution  will  be  obviously  dilTcrcnt  if  the  definition  of  wave  space  is  changed.  For 
instance,  at  any  point  x where  to  and  k arc  interrelated  by  a dispersion  relation 

(I)  li  (k,  A)  (4) 


in  w hich  the  properties  of  the  propagation  medium  arc  characterized  by  the  function  A (x,  t), 
we  can  write  for  the  same  mean  square  surface  deformation 


(5) 
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where  0,  with  k.^)  — (k  cos  0,  k sin  G),  is  the  direction  of  k relative  to,  say,  x and  the 
distribution  9 represents  the  so-called  directional  spectral  density.  and  9 are  connected  by 


(k;  X,  t)  = 


to  til 
k ^ 


9 (o),  0;  X,  t). 


(6) 


Consider  now  a medium  with  nonuniform  still  water  depth  h and  moving  with  velocity  U 
relative  to  the  fixed  (x,  2)  co-ordinate  system.  It  is  assumed  that  the  temporal  and  spatial 
variations  in  h and  U over  any  given  period  and  wavelength,  i.e.  {bhlct)j(i>h  etc., 

are  much  smaller  than  unity,  and  that  the  continuity  equation  for  the  mean  Row  is  satisfied 
(see,  e.g.  Bretherton  and  Garrett,  1969,  p.  554).  We  can  write  (Phillips,  1966,  p.  43) 

(0  = a (k,  3.)  = U . k + 0)'  (7) 

where 

w'  = [gk  tanh  khyi^  (8) 

represents  the  frequency  relative  to  a coordinate  system  moving  with  the  current  U.  Further- 
more, the  equations  for  wave  rays  are  given  as  (see,  e.g.  Bretherton  and  Garrett,  1969; 
Kenyon,  1971) 

dx  cLl  , dk  dil  bX 

df  bk  df  bX  b\ 


Of  the  preceding  expressions,  the  first  describes  a wave  ray,  i.e.,  the  path  traced  out  by  an 
observer  moving  with  the  absolute  group  velocity 


in  which 


^(U.k)  -fC'c 

ck 


C'n  = 


b(o'  k 
bkk 


(10) 

(11) 


denotes  the  group  velocity  relative  to  U,  and  the  second  of  (9)  describes  the  change  in  the 
vector  wave  number,  k,  along  the  ray. 

Under  conservative  conditions  in  which  dissipation,  generation  and  wave-wave  inter- 
actions can  be  neglected,  the  energy  balance  reduces  to  (see,  e.g.,  Willebrand,  1975,  p.  125) 

where  Va  = {bjbx,  blby)  represents  the  horizontal  gradient  operator. 

Clearly  (12)  describes  the  propagation  of  the  quantity  ('.^,  co'),  referred  to  as  wave  action 
spectral  density,  in  a genera'  time  dependent  inhomogeneous  medium,  and  implies  that 

2!.  = ^ = constant  (13) 

0) ' kb  k 10 ' 
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along  a wave  ray.  f hercforc,  thi.s  result  together  with  (9)  and  the  initial  values  of  'y  (or  9), 
X and  k at  one  time  are  suilicient  to  determine  the  speetral  densityy  (k;  x,  /)  or  9 (w,  0;  x.t ). 
An  obvious  alternative  that  avoids  the  excessive  numerical  integration  involved  in  rav 
tracing  in  general  is  to  solve  (12)  by  a finite  dilTerencc  approximation,  this  time,  together 
with  the  irrotationality  condition  and  the  kinematical  conservation  equation  for  the  vector 
wave  number,  k (see,  c.g.,  Phillips,  1966,  p.  4.4). 

In  order  to  explore  more  explicitly  what  is  involved  in  (13)  in  the  following,  we  consider 
the  simple  and  interesting  ca,se  of  steady  state  conservative  random  waves  propagating  from 
a spatially  homogeneous  region  such  as  deep  watei  into  an  inhomogeneous  region  with  a 
nonuniform  depth  profile  h (.v),  and  traversing  a steady  nonuniform  current  field,  L' 
[t/(.v),  K(.v),  lk'(.v,  z)],  as  schematically  illustrated  in  Fig.  I.  Under  these  conditions,  the 


RAY  =.  orthogonal 


Fig.  I.  Schematic  diagram  showing  various  dclinitions  and  (inscrti  ihc  hmil  c.isc  . 12  . t 0 
spatial  inhomogeneity  of  the  problem  is  restricted  only  to  the  v-direction  with 


k sin  0 constant. 


(14) 


and  the  absolute  frequency  to  is  invariant.  The  angle.  (1.  relative  to  the  v-axis  m taken  to  be 
positive  clockwise.  It  is  now  immediate,  from  (7)  and  (13)  with  l’„  ((  ( v).  I ( v),  t)].  th.it 


•ylk) 

(I)  — U/, . k 


(0  f if  9 (dl.  D) 

k (k  0)  — . k 


constant 


(15) 


in  the  jv-direction.  It  is  understood  in  (15)  that  y ? me  independent  of  t.  and  the 
.Y-dependcncy  of  both  quantities  arc  kept  implicit  for  s mplicity.  .AKo.  to  spceifv  various 


1 


I 
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spaec-depciuleiit  quantities  more  explieitly,  yye  will,  from  now  on,  designate  deep  water 
values  that  are  spatially  homogeneous  in  the  absence  of  currents  by  the  subscript  (cjd),  those 
m linite  depth  and  in  the  absence  of  currents  by  the  zero  subscript  and  leave  all  values  in 
the  presence  of  a current  unsubscripted  irrespective  of  any  depth  consideration.  In  this 
manner,  it  is  noted  that  the  simplest  I'orm  of  the  constant  in  (15)  is  or  (Cq 

In  terms  of  these  and  from  (15)  it  follows,  therefore,  that  the  inhomogeneous  densities 
V (k)  and  7 (o),  0)  are  given  by 

V(M  ~ ■;»(k^)  (16) 


9 (<'),  0) 


- ‘I;  '') 

Aoc  (a>/a) 


(<'>.  0,„). 


(17) 


The  general  character  of  the  preceding  results  indicates  that  the  interaction  between  a 
random  wave  field  and  nonuniform  current  depth  elTects  involves  a spatial  transformation 
of  the  spectral  magnitudes,  7^,  and  their  respective  wave  spaces.  However,  before  we 
proceed  to  interpret  this  transformation,  it  is  appropriate  to  discuss  various  kinematical 
and  dynamical  constraints  embedded  in  the  above  equations.  First,  note  from  (16)  and  (17) 
that  the  condition 


(18) 


must  be  salis/icd,  in  particular,  by  Ifie  components  propagating  in  the  direction  of  the  local 
current.  V However,  it  is  also  noted  that  these  components  must  locally  have  |sin  0 | < I . 
T he  latter  constraint  can  be  expressed,  using  (7),  (8)  and  (14),  in  the  form 


(1  > [tanh(kh)|sinO,|]‘« 


(19) 


Recognizing  that  I > tanh  (kh)  |sin0,,|  >0,  it  becomes  eviden'.  that  the  kinematical 
constraint  implied  by  (19),  only,  is  significant.  In  essence,  this  constraint  is  an  extension  of 
that  given  by  Longuet-Higgins  and  Stewart  (1961,  p.  547)  to  the  inclusion  of  the  effect  of 
nonuniform  water  depth.  At  the  lower  limit  where  (19)  becomes  an  equality,  |sin  0|  = I, 
and.  therefore,  we  assume  in  analogy  with  monochromatic  wave  behavior  that  the  associated 
spcctnil  component  is  totally  rellccted  (see,  c.g.,  Longuet-Higgins  and  Stewart,  1961; 
Kenyon,  1971). 

The  quantity  rU/r/i  in  the  denominator  of  (1 7)  represents  the  component  of  the  transport 
velocity,  Q,  along  the  orthogonal  or,  equivalently,  in  the  propagation  direction  (k/A)  of  a 
component  wavelet,  fherefore,  locally,  the  condition 

Q;-''  C\:  > U„.‘‘>0  (20) 

( l<  K k 


I 


I 
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must  be  satislieil.  In  other  words,  the  local  group  velocity,  C'g,  must  be  opposite  in  direction 
and  larger  in  magnitude  relative  to  the  horizontal  current  component,  li//  . (k'A),  in  the 
direction  of  wave  propagation.  In  the  limit  condition  when  C'q  --  — U//  . (k/A:), 
as  schematically  illustrated  in  the  insert  to  Fig.  I,  the  associated  spectral  component  can  no 
longer  propagate  against  the  current  in  that  direction.  Theoretically,  the  local  spectral 
magnitude,  9,  becomes  infinite.  As  in  the  case  of  monochromatic  waves  (see,  e.g.,  Longuet- 
Higgins  and  Stewart,  1961 ).  this  suggests  that  these  spectral  components  will  tend  to  diminish 
or  attenuate  by  uaie  breaking  and,  possibly,  by  a lateral  stretching  in  the  crest  direction 
before  this  point  is  reached.  Note  from  (16)  that  the  corresponding  spectral  magnitude, 
•y,  over  k-space  remains  always  bounded  (i.e.,  < 2).  It  is,  however,  evident  that  (16)  as  well 
as  all  other  results  and  definitions  of  the  preceding  analysis  will  lack  validity  near  the 
critical  point. 

In  summarizing  the  preceding  discussion  now,  we  may  conclude  that  the  spatial  trans- 
formations of  the  spectral  magnitudes,  and  9,^  via  (16)  and  (17)  are  subject  to  the 
reflection  and  breaking  constraints  (19)  and  (20),  respectively.  The  transformations  of  the 
incident  densities,  9^  and  arc,  therefore,  continuous  in  .v  for  all  but  the  attenuated  and 
reflected  components.  The  regions  of  the  incident,  k„  or  w,  0^-space  that  violate  either 
of  these  constraints  at  a point,  .v,  are  cut-off,  or  simply  deleted  beyond  that  point,  assuming 
the  absence  of  any  sort  of  interaction  between  the  attenuated  and/or  reflected  components 
and  tho.se  remaining.  However,  the  original  incident  wave  space  is  to  be  properly  modified, 
as  will  be  illustrated  later  with  examples,  by  the  reflected  components  in  a manner  consistent 
with  the  steady  state  assumption. 

fn  general,  the  spatial  transformation  of  the  magnitude  is  entirely  due  to  the  current 
interaction  in  contrast  with  the  transformation  of  9^  that  involves  the  combined  current- 
depth  effects.  The  spectral  magnitude,  y,  relevant  to  components  with  a locally  opposing 
angle  to  the  current  (i.e.  U„  . k < 0)  are  amplified  in  comparison  with  y^.  and  those 
propagating  with  the  current  (i.e.,  II„  . k ■ 0)  are  suppressed  irrespective  of  any  depth 
consideration.  In  the  case  of  9,  the  current  effect  on  components  locally  opposing  the  current 
is  qualitatively  the  same  as  in  the  case  of  'y,  i.e.,  amplification.  Howevre,  for  components 
propagating  with  the  current,  the  net  effect  of  the  combined  current-depth  interactions 
could  bean  amplification  or  suppression  of  9 relative  to  9^  depending,  from  (17),  on  whether 
the  a.ssociated  to,  0,„  values  satisfy  the  condition 

-AfQ),  (1  M (21) 

rk  \ CO  / 

I 

t 

I or  not,  respectively. 

j The  local  k-space  associated  with  the  transformed  magnitude  y is  distorted  by  the 

t combined  current-depth  effects  in  terms  of  the  magnitude  k as  well  as  the  direction  6 of  the 

i wave  number  vector  k.  On  the  other  hand  the  distortion  in  the  polar  co,  0-space  is  entirely 

I due  to  the  spatial  dependence  of  0.  Since  the  complete  spectral  transformations  require  the 

j mapping  of  the  density  ly  or  9 desirably  in  the  form  of  contours  in  the  local  wave  space,  the 

j invariant  nature  of  the  frequency  (0  makes  the  polar  co,  0-space  particularly  advantageous 

' to  work  with.  Therefore,  the  discussion  will  be  restricted  to  this  space  from  now  on. 
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HORIZONTAL  SHEAR  CURRENT  CASE,  EXAMPLES,  AND  VARIATION 

OF  MEAN  ENERGY 

Considering  the  simplest  case  of  a horizontal  shear  current,  U,,  ==  [0,  K (jc),  0],  and  using 
(7),  (8),  (10),  (11)  and  (14)  result  in  the  following: 


i)''^  r 1 — sin  = gk  tanh  kh 

L d 


represents  the  dispersion  relation,  and 

C'  = C.[l  ^ 

is  the  phase  speed  relative  to  the  current.  Also 
sin  0 = (k^/k)  sin 


y{x)  . „ 
— sm  0 


•]' 


tanh  kh 


where 


= (tanh  <;/i  Fl  — ^ sin  0„ol  sin  0 

L Cgo  -j 

- Fl  — sin  0^1  (tanh  kh)^ 

l_  ^ ao  -I 

] 

» - . fi  i 1 . 

C 2 L sinh  2kh  J 


fii  tor  V '''W  • n 

= n -I  (1  — n)  sm  0 

ck  A:  L 


C',; 


On  substitution  from  these  (13),  (19)  and  (20)  become,  respectively. 


9 (to,  0) 


AC„[l  - sinO^j 

k j^/; 


2 A-  I /;  i (1  — n)  sin  O^^J 


?®  (W>  Ooo) 


(I  - n) 


— r (.v)  sin  0,„,  0„  . 0. 


(22) 


(23) 


(24) 

(25) 

(26) 

(27) 


(28) 


j^l  - sin  0„j  > [(tanh  kh)  sin  0„]'^  0*  > 0 (29) 


(30) 


At  this  point,  it  is  worthwhile  to  discuss  some  e.xamples  illustrating  what  is  involved  in 
the  preceding  transformations  and  constraints.  For  simplicity  in  presentation  the  examples 
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will  be  contined  to  deep  and  shallow  water  conditions  where  various  key  definitions  involved 
in  the  transformation  (28)  become  analytically  tractable.  Furthermore,  it  will  be  convenient 
to  nondimensionalize  various  quantities  in  terms  of  the  value  = max  V (x)  as  follows: 


k* 


gh  (x) 


(31) 


The  explicit  form  of  the  incident  spatially  homogeneous  deep  water  spectral  density  in 
both  examples  is  assumed  to  be 


1 ; 1/2  < w*  < 3,  je^l  < n/2. 
0 ; otherwise 


(32) 


The  incident  (0*,  0^-space  has,  therefore, a semi-annular  shape  as  schematically  illustrated  in 
Figs.  2 and  7 (continuous  line  boundary).  A spectral  density  of  the  form  (32)  by  no  means 


Fto.  2.  Example  I : Current  interactions  in  deep  water.  Definition  sketch  showing  the  qualita- 
tive effect  of  the  shearing  current  on  wave  orlhogonals  for  various  angles  of  entry  9qo. 


characterizes  any  oceanic  situation.  This  artificial  form  is  of  convenience  here  in  that  the 
deformation  of  the  incident  w,  0^-space  and  contours  of  density  amplification  or  suppression 
can  be  followed  easily.  Also,  because  of  the  unity  magnitude,  the  transformation  of  a more 
realistic  spectral  density  is  readily  obtained  with  a simple  multiplication  of  amplitudes  over 


regions  of  co,  0„-space  coincident  with  the  one  chosen  here. 

Example  I : Current  interactions  in  deep  water 

In  deep  water  (24),  (29),  (30)  and  (28)  become,  respectively, 

sin  0 = fl  — CO*  sin  0„"|  sin  0,^  (33) 

L T 

(0*  sin  0„  < (1  - sin"2  0,„),  0„  > 0 (34) 

^ m 

V (x) 

- CO*  sin  0„  < 1,  0„  < 0 (35) 

9*  (co*,  ((0*,  0„^)  = A (36) 

where  the  spectral  amplification  factor,  A is  in  the  form 

^ = fl  - (0*  sin  0„]T  1 4 CO*  sin  0„1  ''  (37) 

L J L m A 


The  shear  current  field,  V (x),  shown  schematically  in  Fig.  2,  consists  of  regions  of  mono- 
tonic increase  and  decrease.  It  is  evident  from  (37)  and  (33)  tha 

> 1 (e„  < 0) 

= 1 (e„  = 0)  (38) 

< 1 (0«  > 0) 

and  the  contours  of  density  (p*  = A = constant,  corresponding  to  various  fixed  values 
of  the  product  co*  sin  0„,  are  straight  lines  parallel  to  the  direction  of  normal  incidence 
(0(0  =0  =0°)  over  both  o)*,  0„-  and  to*,  0-space.  The  breaking  constraint  (35)  indicates 
the  regions  of  the  incident  co*,  0„-space  to  be  deleted  in  a progressive  manner  dependent 
on  the  ratio  V (x)/F„  and  whether  dK/dx  > 0 or  dK/d.x  < 0.  For  instance,  referring  to 
Fig.  2,  the  spectral  components  in  the  region  ABL  of  the  incident  co*,  0„j,-space  are  entirely 
eliminated  by  the  time  waves  propagate  to  the  section  (2-2)  where  V (.x)  = 0.5  V,,.  Simi- 
larly, the  components  in  the  region  BCKL  gradually  dissipate  as  waves  advance  from  (2-2) 
to  (3-3).  However,  beyond  (3-3)  where  dV/dx  < 0,  wave  breaking  has  no  influence  on  the 
spectral  components  to  the  right  of  the  line  CK.  For  0„  > 0,  the  reflection  constraint  (34) 
suggests  that,  as  waves  propagate  from  (l-l)  to  (2-2),  the  part  IFGH  of  the  incident  co*. 
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*'oc'SPdce  be  eliminated  at  (2-2).  Similarly,  from  (2  2)  to  (3  3),  the  components  in  the 
region  FiFIJ  are  reflected  leaving  only  the  part  CDf£JK  of  the  incident  wave  space  beyond 
(3  -3)  where  breaking  and  reflection  become  immatcriaf.  However,  for  a steady  state,  the 
retlected  spectral  components  must  be  properly  accounted  for  since  they  propagate  back  in 
the  direction  (rt  0)  locally  and  modify  the  incident  wave  space.  Therefore,  tlie  deep  water 
density  7sc*(  0 ^nd  the  associated  lo*,  0^^, -space  must  be  mouilied  as  shown  in  Fig.  2 

by  dash  lines  with  the  region  Ge'J'H  representing  the  mirror  image  of  the  part  GEJH  that 
has  been  progressively  reflected  back  as  waves  propagated  from  (l-l)  to  (3-3). 

Figures  3 5 correspond  to  a simple  mapping  of  the  transformed  densities  over  their 
local  (I)*,  0-spaces  at  any  point  along  the  lines  (2-2),  (3-3)  and  (4  4)  of  Fig.  2,  respectively. 
By  the  lettering  in  Fig.  2 and  the  corresponding  ones  in  these  figures,  the  local  distortion 
of  the  various  regions  of  the  incident  wave  space,  and  the  amplification  or  suppression  of  the 
spectral  magnitudes  are  easily  observed.  It  is  noted  in  each  case  that  the  spectral  amplitude 
of  components  with  normal  incidence  to  the  current  (0,^  - 0 ) traverse  the  current 
unaffected,  whereas  those  with  opposing  (following)  angles  of  entry  are  amplified  (sup- 
pressed). In  Fig.  3,  the  region  FE'J  '1  represents  those  components  reflected  back  in  between 
(2-2)  and  (3-3).  Finally,  Fig.  6 illustrates  the  spectral  density  and  the  associated  wave  space 
at  any  point  on  and  beyond  the  line  (5-5)  where  V (x)  — 0.  The  spectral  magnitudes 
(=  I)  and  the  wave  space  in  this  figure  are,  therefore,  exactly  the  same  as  the  CDEJK  part 
of  Fig.  2. 

0° 


L K 


1 


Fio.  4.  Example  I : Current  interactions  in  deep  water.  Contours  of  <p*  — .4  " const,  and  the 
associated  to*,  9-space  at  any  point  along  section  3-3. 


I 

I 


Fio.  5.  Example  I : Current  intereactions  in  deep  water.  Contours  of  tp*  - A — const,  and  the 
associated  o>*,  6-space  at  any  point  along  section  4 4. 
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Fio  6 Example  1:  Current  intertetions  in  deep  water.  The  dir<xiional  spectral 
, ■ = 1^,  vithin  the  dashed-line  •-  'undary)  and  the  as_s<viated  u.  . 0-space 

nt  ;inv  noint  on  aiid  beyond  section  5 5. 


density 
0)*,  Oso) 


Example  2:  Current  interactions  in  shallow  ii-aler  , • - 

By  virtue  of  the  delinitions  (31),  we  may  rewrite  the  dispersion  relation  ( 


22)  as: 


I 1 


tanh  k*li*  == 


Hy) 

k* 


to*  sin  0 


-T 


(39) 


Under  the  shallow  water  condition 


h*  {(o*f  r 1 - w* 

L ^ m 

sin  0OO  j ’i/'O 

(40) 

we  find  [see,  e.g.,  Longuet-Higgins,  1956] 

tanh  k*li*  = to*  (h*) 

1 — (0*  sin  0„1 

v,„  J 

(41) 

and 

^ (tanh  k*h*)  ^ 2 w*  (/i*)''M 

tk*  L 

1 — w*  sin  0„ol  . 

y,n  J 

(42) 

It  is  immediate  from  these  approximations  that 

C — C (j  — 

Vn,  (/'*)*'“ 

(43) 

and  (24),  (29),  (30)  and  (28)  become,  respectively. 
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w*  sin  0^  -r  [1  4-  (//♦)> > 0 (45) 

^ m 

- ^ w*  sin  ~ 00,  0„  < 0 (46) 

y,„  I - n 

A=<f*  (o)*,  0)/9„*  (to*,  0^)  =.  I I - to*  sin  (2//*  ((o)*‘‘]->.  (47) 

As  a specific  example  here,  consider  a current  of  the  form  schematically  illustrated  in 
Fig.  7,  together  with  the  shallow  water  depth  profile.  We  see  that  the  breaking  constraint 


Fig.  7.  Example  2:  Current  interactions  in  shallow  water,  netinillon  sketch  in  plan  view  (lop) 

and  profile. 


(46)  is  redundant  in  this  case.  That  is  to  say,  no  component  with  finite  to*  in  the  local  to*, 
Og-space  just  before  the  current  (see  Figs.  7 and  8)  will  ever  attain  a breaking  condition  on 
account  of  excessive  depth  refraction  which  distorts  and  focuses  the  wave  space  in  the 
direction  of  normal  incidence  to  the  current.  Note,  from  (15)  and  (17)  with  IJ|/  = 0,  that 
the  spatial  modification  of  spectral  amplitudes  due  to  depth  effects  is  described  by 

/<='''*=  (C'g/^)oo 
9 (Q/^) 


i 
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0° 


Fig.  8.  Exaniple2:  Current  interactions  in  shallow  water.  Diagram  of  contours  of  ((>„*  A 

const,  and  the  associated  to*,  0„-space  at  any  point  along  the  line  (1-1),  showing  the  dramatic 
effect  of  pure  depth  refraction  on  the  incident  (ABCDEF)  spectral  components  as  well  as  on 
those  progressively  reflected  back  (DF'E)  by  the  shear  current  in  between  (l-l)  and  (2-2)  of 

Fig.  7 (A*  0.004), 


Evidently  (48)  represents  a slight  suppression  (i.e.,  A < 1)  at  first  as  waves  feel  the  bottom, 
then  becomes  a monotonically  increasing  amplification  (A  > 1),  which  is  of  the  form  (47) 
with  y (x)  = 0 in  shallow  water.  However,  this  characterization  holds  in  general  depending 
on  the  local  water  depth  only  rather  than  on  the  nature  of  isobaths  or  any  directional 
consideration. 

The  lowest  upper  bound  of  the  reflection  constraint  (45)  (w*  sin  0„,  =:  0.9572)  is  realized 
at  points  along  the  line  (2-2)  beyond  which  no  reflection  takes  place.  This  bound  is 
schematically  illustrated  with  the  dash  line  DF  in  Fig,  7.  Therefore,  as  waves  reach  (2-2), 
the  components  in  the  DFE  region  of  u)*,  0„-space  are  progressively  reflected,  propagating 
back  and  locally  modifying  the  incident  wave  field  as  in  the  previous  example.  Correspond- 
ingly, the  steady  state  deep  water  wave  space  would  be  as  shown  in  the  same  figure,  with  the 
region  DF'E  bounded  by  dash  lines  representing  the  mirror  image  of  DFE.  Presented  in  the 
subsequent  Figs.  8-1 1 are  the  local  spectral  contours  and  wave  spaces  corresponding  to  any 
point  (l-l),  (2-2),  (3-3)  and  (4-4)  of  Fig.  7,  respectively.  In  particular.  Fig.  8 illustrates  the 
depth  affects  discussed  in  the  preceding  on  the  incident  and  reflected  spectral  components 
just  before  entry  to  the  shear  current.  Figures  9 and  10  present  the  combined  current-depth 
effects.  Qualitatively,  components  with  an  opposing  angle  of  entry  (0,^,  < 0)  are  further 
amplified  and  focused  in  the  direction  of  normal  incidence.  Obviously,  for  components 
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Fio.  9.  Example  2:  Current  interactions  in  shallow  water.  Contours  of  tp*  = /(  = const,  and 
the  associated  u)*,  6-space  at  any  point  on  the  line  (2-2),  showing  the  combined  effects  of  depth 
and  current  interactions  (A*  = 0.002). 


0" 


Flo.  10.  Example  2:  Current  interactions  in  shallow  water.  Contours  of  tp*  .-1  = const, 
and  the  associated  O)*,  6-space  at  any  point  on  the  line  (3-3),  showing  again  the  combined  effects 
of  depth  and  current  interactions  (A*  0.001). 
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Fig.  II.  Example  2:  Current  interactions  in  shallow  water.  Contours  of  <p„*  = .-1  = const, 
and  the  associated  to*.  0j-space  at  any  point  on  the  line  (4-4)  and,  beyond  showing  the  effect  of 
pure  depth  refraction  right  after  waves  cross  the  current  (A*  = 0.001). 

following  the  current  (0„  > 0),  current  and  depth  interactions  have  an  opposing  influence 
and  their  net  effect  is  in  general  characterized  by  (47).  Therefore,  local  spectral  amplitudes 
for  which  (o*,  0„  values  satisfy  (violate)  the  condition  2h  ((«*)''“  < [1  — (f^/T'/n)  w*  sin 
remain  amplified  (suppressed)  relative  to  the  corresponding  amplitudes  in  deep  water  as 
illustrated  in  Fi"  > to  the  left  (right)  of  unity  contour  line.  Components  dominated  over  by 
the  current  (i.e.,  A < 1}  are  in  effect  swept  along  with  the  current  as  they  propagate  from 
(1-1)  to  (2-2).  Between  (2-2)  and  (3-3)  where  the  current  is  uniform,  only  depth  effects  are 
operative  in  amplifying  and  focusing  spectral  components  further  relative  to  those  at 
(2-2).  This  is  illustrated  for  any  point  along  (3-3)  in  Fig.  10.  Finally,  beyond  (3-3)  where  the 
depth  is  uniform,  only  current  effects  are  of  concern,  interacting  with  the  wave  field  in  a 
reverse  and  diminishing  manner,  exactly  as  in  the  region  from  (3-3)  to  (5  -5)  of  the  previous 
example.  Therefore,  spectral  contours  and  wave  space  would  be  as  shown  in  Fig.  1 1 at  any 
point  on  and  beyond  (4-4). 

Comments  on  mean  energy.  Mean  total  wave  energy  per  unit  horizontal  area  is  defined  by 
(Phillips,  1966,  p.  27) 


£ = pg  < T)-  > = pg  I 9 (w,  0)  u)  du)  dO.  (49) 

C0,O 

In  order  to  comment  briefly  on  the  spatial  variation  of  this  quantity  for  the  particular  case 
of  interest,  here,  it  will  be  expedient  to  write 


cO 

9 (o),  0)  to  dw  dO  = 9 (o),  0)  u)  dw  d6„ 


where,  by  virtue  of  (24), 


^0 


k cos  0 


sin  0„  ck  \ 

BQj 


(50) 


(51) 


and,  from  (25), 


I 

I 
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i k 


k k„  k'(X)  cos  0„ 

/)  (o  1^1  — ^ sin  0„  I 


On  substitution  from  (50),  (51)  and  (52),  (49)  can  be  written  as 

i-'  = Pg  f 9x  0®)  w dm  de„ 

J n sin  20 


0).  0 _ e K 


(52) 


(53) 


where  R represents  the  region  of  the  incident  to,  0„„-space  that  is  not  deleted  by  wave 
breaking  and/or  reflection  as  waves  propagate  to  the  point  of  interest.  This  expression  is  a 
generalization  of  various  well-known  results  associa'.ed  with  monochromatic  waves  and, 
of  course,  with  a narrow-band  random  wave  field.  For  instance,  in  the  case  of  mono- 
chromatic waves  with  the  incident  height,  //„,  frequency  w,  direction  0^,  and  energy 
---=  (l/8)pg  one  can  show  from  (53)  that 


(n  sin  20)„  ^ 
n sm  2U 


(54) 


This  obviously  reduces  to  the  familiar  form 


^ (Cc  cos  0)^,  ^ 
(Co  cos  0)o 


(55) 


in  the  absence  of  any  currents  (see,  e.g.,  Phillips,  1966,  p.  53),  and  to 


^ sin  20a,  ^ 
sin  20 


CO 


(56) 


in  deep  water  with  the  presence  of  currents  (Longuet-Higgins  and  Stewart,  1961,  p.  547). 


CONCLUDING  REMARKS 

A mean-square  spectral  distribution  constitutes  a concise  characterization  for  incoherent 
random  gravity  waves.  In  the  preceding,  the  spatial  transformation  of  this  distribution  by 
refraction  due  to  currents  and  bottom  topography  has  been  presented  in  closed  form  and  with 
explicit  examples  in  the  particular  case  of  one-dimensional  variations  in  currents  and  topo- 
graphy. Clearly,  the  omission  of  various  other  realistic  effects  such  as  nonlinear  mechanisms, 
frictional  dissipation,  generation,  etc.,  requires  a tentative  judgment  on  the  results  here 
as  in  the  case  of  linear  monochromatic  waves.  It  is  seen,  however,  that  currents  and  under- 
water topography  can  have  a dramatic  influence  on  the  spectral  characteristics  of  random 
waves  Shear  currents,  in  effect,  act  as  a filter,  dissipating  and/or  reflecting  certain 
o-mponents  while  transmitting  others  with  a substantial  amplification  or  reduction  in 
•^pliiudc  These  interactions  are,  then,  modified  further  and  become  more  complex  in 
)'  'wer  depths  with  depth  refraction  and  shoaling.  Various  aspects  of  these  interactions. 
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particularly  kinenialics,  arc  well  known  in  the  case  ol  deterministic  monochromatic  waves. 
The  primary  motivation  in  the  extension  of  these  concepts  to  incoherent  random  waves 
lies  in  the  more  general  character  of  results  in  terms  of-a  spatially  inhomogeneous  spectral 
distribution  and,  therefore,  in  the  fact  that  it  supplements  the  statistical  approach  to  the 
description  of  the  real  ocean  surface.  This  approach  is  presently  considered  by  many  as  a 
more  fruitful  one  than  a deterministic  theory  based  on  a monochromatic  wave  train. 

Financial  support  for  this  work  was  provided  by  the  Geography  Programs  of  the  Office 
of  Naval  Research  under  contract  No.  NOO 14-76-0-0412  with  the  University  of  Delaware.  The  authors 
gratefully  acknowledge  the  helpful  comments  by  the  referees  and  R.  G.  ITean. 
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model  taking  into  account  the  varying  water  depth  in  shallow  water.  The  third 
part  presents  a method  of  calculating  the  two  dimensional  random  wave  field 
near  shore  in  terms  of  directional  spectrum,  taking  into  account  the  wave- 
current  interactions. 
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